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A Correlated Study of Optical and X-ray Afterglows of GRBs 
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ABSTRACT 

We study an extensive sample of 87 GRBs for which there are well sampled 
and simultaneous optical and X-ray light-curves. We extract the cleanest possible 
signal of the afterglow component, and compare the temporal behaviors of the X- 
ray light-curve, observed by Swift XRT, and optical data, observed by UVOT and 
ground-based telescopes for each individual burst. Overall we find 62% GRBs 

department of Physics, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden; fryde@kth.se 
2 The Oskar Klein Centre for Cosmoparticle Physics 

^Department of Physics, Stockholm University, SE-106 91 Stockholm, Sweden; liang.li@fysik.su.se 
4 Erasmus Mundus Joint Doctorate in Relativistic Astrophysics 

5 Chinese Center for Antarctic Astronomy, Chinese Academy of Sciences, Nanjing 210008, China 

6 Joint Center for Particle Nuclear Physics and Cosmology of Purple Mountain Observatory and Nanjing 
University, Nanjing 210008, China 

7 Key laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, Chinese Academy 
of Sciences, Nanjing 210008, China. 

department of Astronomy, Nanjing University, Nanjing 210093, Jiangsu, China 

9 Key laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, 
Nanjing 210093, China 

10 GXU-NAOC Center for Astrophysics and Space Sciences, Department of Physics, Guangxi University, 
Nanning 530004,China 

11 Department of Physics and Astronomy, University of Nevada, Las Vegas, NV 89154, USA 

12 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China 

13 Center for Space Plasma and Aeronomic Research (CSPAR), University of Alabama in Huntsville, 
Huntsville, AL 35899, USA 

14 Dip. di Fisica and ICRA, Sapienza Universit di Roma, Piazzale Aldo Moro 5, 1-00185 Rome, Italy. 
15 ICRANet, Piazza della Repubblica 10, 1-65122 Pescara, Italy. 

16 National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, China 



2 


that are consistent with the standard afterglow model. When more advanced 
modeling is invoked, up to 91% of the bursts in our sample may be consistent 
with the external shock model. A large fraction of these bursts are consistent 
with occurring in a constant interstellar density medium (ISM) (61%) while only 
39% of them occur in a wind-like medium. Only 9 cases have afterglow light- 
curves that exactly match the standard fireball model prediction, having a single 
power law decay in both energy bands which are observed during their entire 
duration. In particular, for the bursts with chromatic behavior additional model 
assumptions must be made over limited segments of the light-curves in order for 
these bursts to fully agree with the external shock model. Interestingly, for 54% 
of the X-ray and 40% of the optical band observations the end of the shallow 
decay (t °' 5 ) period coincides with the jet break (t p ) time, causing an abrupt 
change in decay slope. The fraction of the burst that consistent with the external 
shock model is independent of the observational epochs in the rest frame of GRBs. 
Moreover, no cases can be explained by the cooling frequency crossing the X-ray 
or optical band. 

Subject headings: methods: statistical — reference systems — X-rays: ISM 


1. Introduction 

While the prompt, gamma-ray emission of gamma-ray bursts (GRBs) generally is in¬ 
terpreted as being due to energy dissipation processes internally in the relativistic outflow 
(e.g. Usov 1992; Rees & Meszaros 1994; Thompson 1994; Kobayashi et al. 1997; Daigne & 
Mochkovitch 1998; Drenkhahn & Spruit 2004; Giannios & Spruit 2006; Zhang & Pe’er 2009; 
Zhang & Yan 2011), the afterglow emission is interpreted as synchrotron emission generated 
in a forward shock as the blastwave interacts with the circum-burst medium (e.g., Meszaros 
& Rees 1997; Sari et al. 1998; Panaitescu et al. 1998; Huang et al. 2000; Panaitescu & 
Kumar 2001; see Zhang & Meszaros 2004 for review). Sari et al. (1998) calculated the broad¬ 
band spectrum and corresponding light-curves and reported quantitative relationships (the 
standard afterglow model) with a constant interstellar medium (ISM), which can generally 
describe the X-ray and optical afterglow emission to be from the same synchrotron emission 
process. Prior to the launch of the Swift satellite observational studies largely confirmed the 
standard model. By comparing the light-curve decay at different wavelengths, Wijers, Rees, 
& Meszaros (1997) showed that the standard model fitted the afterglow data for GRB970228 
and 970402. Likewise, Panaitescu et al. (2005a, b) concluded that the X-ray and optical 
afterglow emission of a sample of 10 bursts were consistent with a single component origin. 
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Furthermore, Nardini et al. (2006) studied the spectral energy distribution over the optical 
and X-ray bands of 24 bursts and concluded that the majority of the bursts were consistent 
with a single emission component. Furthermore, Li et ah (2015 in preparation) recently 
studied the evolution of color indices for a large GRB sample which have well multi-band, 
optical afterglow observations, and also find that most of these GRBs can be explained by 
the standard synchrotron afterglow model. 

Although the standard afterglow model is successful in explaining the overall features 
of the observed afterglow behavior, the increasing amount of observational data, as well as 
its gain in detail, clearly call for additional assumptions to be made, e.g., on the ambient 
density profile, prolonged energy injection, jet, accretion around central engine, and viewing 
angle effects (e.g. Meszaros et al. 1998; Dai & Lu 1998, 2002; Rees & Meszaros 1998; Rhoads 
1999; Sari et al. 1999; Chevalier & Li 1999, 2000; Huang et al. 2000; Zhang & Meszaros 
2001; Pisani et al. 2013; Zhang et al. 2014; Ryan et al. 2014; RufRni et al. 2014). Findings 
based on these assumptions indicate that the afterglow emission is more complicated than 
that predicted by the simplest standard afterglow model, and further developments studying 
“post-standard effects” were made by various authors. For example, Chevalier & Li (1999) 
took a wind-like medium into account, that is, a circum-burst medium density that decreases 
as r -2 , where r is the distance from the center, due to the wind from the progenitor. This 
gave rise to new, quantitative “closure relations” of the external shock model. Detailed 
modeling of early and late afterglow properties in a wind medium has been carried out by 
Panaitescu & Kumar (2000), Wu et al. (2003), Kobayashi & Zhang (2003), and Kobayashi 
et al. (2004). The energy injection models have been studied by Dai & Lu (1998), Rees & 
Meszaros (1998), Sari & Meszaros (2000), and Zhang & Meszaros (2001). The jet models 
are studied in detail by Rhoads (1999), Sari et al. (1999), Huang et al. (2000) and Wu et 
al. (2004), and more recently, simulations have been carried out by Eerten & MacFadyen 
(2013) and Duffell & MacFadyen (2014). Zhang et al. (2006) derived the “closure relations” 
of relevant external shock models and collected them in their Table 2. A more complete 
review of closure relations can be found in Gao et al. (2013). 

The Swift satellite (Gehrels et al. 2004) has provided an unprecedented sample of 
long time series of broad energy band observations of bursts. Simultaneous observations 
with the Swift instruments BAT, XRT, UVOT, and ground-based optical telescopes have 
revolutionized our understanding but have also revealed some critical problems with the 
standard model (e.g., Zhang 2007, 2011). In contrast to the expected results from the 
standard model, many bursts showed a chromatic behavior between that of the X-ray and 
optical band, that is, showing different behaviors at different frequencies. For instance, 
Panaitescu et al. (2006) studied 6 Swift bursts that exhibited distinctly different behaviors 
of the X-ray and optical light-curves, prompting them to argue for a different origin of the 
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emission in the two bands, at least for the bursts in the study. Moreover, Oates et ah 
(2009) studied a sample of 26 Swift bursts and showed that the XRT X-ray and UVOT 
optical light-curves have similar statistical properties. However, they found that they are 
remarkably different during the first 500s after the BAT trigger. The late time light-curves 
have a better resemblcnce. Zhang et al. (2007) and Liang et al. (2007, 2008, 2009) in a series 
of papers, who studied X-ray afterglow and their consistency with the external shock model, 
also suggest that the X-ray and the optical band have a different origin within the emission. 
Systematic studies of Swift X-ray afterglows and their consistency with external shock models 
have been also carried out by O’Brien et ah (2006), Willingale et ah (2007), Evans et ah 
(2007, 2009) and Zaninoni et ah (2013). In (Rufhni et al 2014), GRBs originating from 
binary system is investigated, one family of GRBs is classified as that a standard black hole 
is formed by a neutron star accreting to the critical mass of gravitational collapse, this family 
holds typical scaling laws and overlapping features in the X-ray afterglow (Pisani et al 2013). 

In this paper, we analyze a sample of 87 bursts, focussing on the clean afterglow emis¬ 
sion after omitting the components which are related to internal shock origin, afterglow 
reverse shock emission, and the late-time associated supernova component. These bursts 
are observed by XRT and simultaneously by UVOT and/or ground-based telescopes and 
are able to make a large statistical study. We examine whether, and how often, the X-ray 
and optical afterglow light-curves are consistent with the external shock model. This paper 
is organised as follows. We define the data sample, the fitting results are presented and 
analysed in section 2. We carry out a detailed statistical analysis in section 3. In section 4, 
we perform several tests to the external shock models for all our afterglow samples. A more 
detailed investigation regarding how well the models can interpret the data can be found in 
Wang et al. (2015). We discuss the properties of afterglow break behaviors in section 5, and 
Conclusions and Discussion in section 6. 


2. Sample Selection and Light-Curve Fitting 

We compiled a data base of 270 GRBs with optical afterglow detections made from 
February 1997 until December 2013 which have been presented in GCN Circulars and in 
the literature. We included observations made by Swift/UVOT, as well as other ground 
based telescopes, such as the Gamma-Ray burst Optical/NIR Detector (GROND; Greiner 
et al. 2007, 2008), ROTSE-lll (Akerlof et al. 2003), TAROT (Klotz et al. 2009), RAPTOR 
(Vestrand et al. 2004), and REM (Zerbi et al. 2001). For all these bursts we also compile 
the full sample of XRT data from the Swift data archive. 

In order to study the afterglow emission expected from the external shock model, we 



need to define a subsample of these bursts, since multiple emission components can exist. 
For instance, Zhang et al. (2006) summarized the X-ray afterglow behavior with a synthetic 
cartoon of an X-ray light-curve based on the observational data from the Swift archive, 
showing that X-ray afterglow included five emission components that have distinctly different 
physical origins or processes. Likewise, Li et al. (2012) summarized the optical afterglow 
behavior with a synthetic cartoon optical light-curves, based on an extensive analysis of the 
optical data statistics, morphology and theory. In contrast to the X-ray afterglow, they 
identified eight clear emission components in the optical afterglow light-curves, which may 
have different physical origins or processes (Li et al. 2012). To extract the cleanest possible 
signal of the afterglow component we therefore apply the following 4 rules: 

Rule 1: We exclude the following components (i) the steep decay in early X-ray after¬ 
glow light-curves, that is involved in the tail of prompt emission related to internal energy 
dissipation (Zhang et al 2006) (ii) the flare behavior that may be related to the late-time 
activity of the central engine (Burrows et al. 2005; Fan & Wei 2005; Zhang et al. 2006; Dai 
et al. 2006; Proga & Zhang. 2006; Perna et al. 2006) (iii) Using the relation a > 2 + /?, we 
exclude the ’internal plateau’ samples (Troja et al. 2007; Liang et al. 2007) (iv) the early 
optical reverse shock emission (v) the GRBs-SNe bumps that are typically detected in the 
optical afterglow light-curves at late time. 

Rule 2: We restrict our sample to include only bursts for which there are good, simul¬ 
taneous observations in both the X-ray and optical bands, since we want to compare the 
optical and X-ray light-curve behaviors. 

Rule 3: All the optical light-curves correspond to R-band observations. For bursts which 
only have, or have better, observations in other bands, we corrected the observations into 
the R-band, by assuming a spectrum F v = F u0 u~^°, /3 0 is the optical spectral index. This is 
particularly the case when a longer time series exists in other bands. 

Rule 4 : We fit the broadest possible time interval with either a single power-law or a 
broken power-law. Such a fit will ignore wiggles on short time scales. For instance, GRB 
051109A is fitted with a broken power-law, which ignores the flare behavior (maybe a late 
re-brightening bump behavior) that is embedded in the light-curve in the narrow time range 
[10 4 s, 10 5 s], 

These rules give us a total of 87 bursts which we classify as four samples depending 
on their light-curves in the X-ray and optical bands. We denote these samples as Group I, 
Group II, Group III, and Group IV which are defined in the following way: 

• Group I: All the bursts for which there is a clear light-curve break detected in both the 
X-ray and optical afterglows, as shown in the left-hand panel in the cartoon of light 
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curves in Figure |TJ Of 87 GRBs, our full sample consists of 24 such bursts, and the 
broken power-law fits to the light-curves are shown in Figure [3] 

• Group II: All bursts for which there is a break detected only in the X-rays, while the 
optical band generally shows a single power-law decay. Note that a few bursts have a 
detected afterglow onset which is followed by the characteristic decay of the sample. 
A cartoon of a typical light-curve in this sample is shown in the centre-left panel in 
Figure [0 Of 87 GRBs, our full sample consists of 24 such bursts, and the fits to the 
light-curves are shown in Figure [4] We note that Group II GRBs do not always point 
towards a chromatic behavior. This is because in many GRBs there is no optical data 
around the X-ray break time, so that the existence of an X-ray break around the same 
time is not ruled out. This issue is fully addressed in Wang et ah (2015). 

• Group III: All bursts for which the break is detected only in the optical band, while the 
X-ray light-curve in general is described by a single power-law. Note that a few bursts 
have a detected afterglow onset. A cartoon example is shown in the centre-right panel 
in Figure [D Of 87 GRBs, This sample comprises 9 bursts and fitted afterglows are 
depicted in Figure [5] Again some GRBs in this group do not have X-ray data around 
the optical break time, so that they are not necessarily chromatic (Wang et al. 2015). 

• Group IV: All bursts which do not have any obvious breaks, either in the X-ray or 
in the optical bands. They are fitted with single power-law functions. A few of these 
bursts have afterglow onsets. Of 87 GRBs, this sample includes 30 bursts and the best 
fitted light-curves are shown in Figure [6] 

One example of the light-curve fittings is given in Figure [2] 

These light-curves were fitted by either a single power-law (SPh) 

F = F 0 t~ a , (1) 

or a smoothly broken power-law (BKPL) 

F = F 0 [(t/t h ) a '“ + (t/t h ) a ^]- 1/u} - (2) 

Here a, oq, «2 are the temporal slopes, t b is the break time, F), = F 0 2~ 1 F is the flux of 
the break time, u describes the sharpness of a break; the smaller the oj, the smoother the 
break. We often fixed this parameter to 3 (for a few cases which have a smoother break we 
fixed it to 1 which improves the fit). Equation ([2]) is used both for light-curves which have 
a break and light-curves which exhibit the afterglow onset. These functions are fitted to the 
data with a non-linear, least-square fitting using the Levenberg-Marquardt algorithm. We 
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use an IDL routing called mphtfun.prcpj (Markwardt 2009). Errors in the parameters are 
calculated either through error propagation or through the bootstrap method, for which we 
produce 1000 realisations which give a Gaussian distribution, whose standard deviation is 
an estimate of the parameter error. 

The fitting results of the X-ray afterglow, including the break time (T Xi b), the power- 
law decay slopes before (a x ,i_) an d after (a x>2 ) the breako, and time intervals during which 
model fittings are performeqj, [T Xjl , T Xj2 ], are summarized in Table 1. Based on the fitting 
parameters, we also calculated the fluence of the X-ray band (S'*), which is the integral of 
the observed flux over the time internals. We also calculate the properties in the rest frame 
for bursts that have a measured red shift (80 of 87 bursts), including the isotropic luminosity 
at the break (L Xib ) and the total isotropic energy release in the X-ray band (A Xjiso ). The 
photon spectral index before (Tvi) and after (T Xj2 ) the break time, are also given and were 
retrieved from the Swift websitcl- Most of the spectral indices were derived using the PC 
mode. Our results are reported in Table 1. 

Similarly for the X-rays, in Table 2, we report the fitting results of the optical afterglow 
over the time interval [T 0)1 , T 0)2 ], the break time, T 0>b , the power-law slope before (after) the 
break, a Q) \ (a 0 , 2 ), the fluence in the optical band, S Q , the break luminosity, A 0 ,b, and the 
isotropic energy release in the optical band, E^ )iso . We also include in the table the optical 
afterglow spectral index, (3 0 , and the redshift, z, found in the literature. We converted the 
observed magnitudes to energy fluxes and calculated the magnitude of the A'-correction by 
k = —2.5 (/3 0 — 1) log(l + z) while the Galactic extinction correlation was calculated by a 
reddening map presented in Schlegel et al. (1998). Any flux contribution in the very late 
epoch (~ 10 6 s after the GRB trigger) that possibly comes from the host galaxy is also 
subtracted. 

A comparasion of the relation properties between the X-ray and the optical afterglow is 
summarized in Table 3. We compare the X-ray and optical data set and, in particular, the 
fit results of the light-curves. We present the ratio of break times (tA^-). Next, we calculate 

-*o,b 7 


1 http: //purl.com / net/rnpfit 

2 Notice that for the light-curves with a single power-law decay we use subscript ”2”. For example, the 
X-ray and optical temporal indices are denoted as a Xj 2 and a 0 , 2 - Similarly, the X-ray photon spectral indices 
in Group III and IV samples are denoted as r x 2 , since the X-ray light curves in these samples are single 
power laws. 

3 The time intervals are selected according to the four rules listed above. We use the subscripts ”x” and 
” 0 ” for parameters in the X-ray and optical bands, respectively. 

4 http : // www.swift.ac.uk / xrt_curves/allcurves.php. 




the differences in power-law slopes in different temporal segments and/or different energy 
bands as well as differences in the spectral slopes in the two bands. We define 
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and present their values in Table 3. We also calculate the X-ray-to-optical ratio of the break 
luminosity (yyy); isotropic energy release ( ^*’[ so ), as well as the flux at 1 hour (^). 


3. Comparison of the Statistical Properties of the X-ray and Optical Bands 

Totally 833 Swift-GRBs were detected by the Swift Satellite until the end of 2013 (from 
December 18, 2004 to December 31, 2013), of which only one-tenth (87/833) had good, 
simultaneously observed X-ray and optical afterglow data. The fraction is far smaller than 
of the totally 691/833 Swift-GRBs with good XRT data observations. 

The distributions of statistical properties comparing the X-ray with the optical band 
are shown in Figure [7] We compare the distribution of decay slope of the X-ray and the 
optical bands in Figure [7] (a). We gather all information on X-ray and optical decay slopes. 
Totally we had 135 decay slope indices for the X-ray (ct x ) and 120 for the optical bands (ct 0 jE 
The distribution of ct x ranges within [-0.01, 2.89] and of ct Q within [-0.08, 2.59]. The X-ray 
emission usually decays faster than the optical emissiorfl. This is consistent with previous 
findings (e.g. Oates et al. 2011; Zaninoni et al. 2013), and is consistent with theoretical 
predictions of the afterglow model (e.g. Sari et al. 1998; Gao et al. 2013). In Figure [7] 
(b), (c), (d), (e), we show the distributions of the break times and the break luminosity for 


including the pre- and post-break segment of all Group samples. X-ray statistics including temporal 
slopes: 48 ct x ,i (Group I+Group II) and 87 a x ^- (Group I+Group II+Group III+Group IV). Optical 
statistics including: 37 ct 0 ,i (Group I+Group III) and 87 ct 0 , 2 : (Group I+Group II+Group III+Group IV). 

6 As summarized in Table 3, overall the Act value of 75 out of 111 (68%) GRBs are greater than 0 (13 of 
24 cases for Acti and 62 of 87 cases for Aa^), which indicates that the X-rays usually decay faster than the 
optical. 
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the bursts which have a break behavior in their afterglow light-curves L b; X-ray: Group 
I and Group II; Optical: Group I and Group III), of the fluences (S X ,S Q ), and finally the 
distribution of the isotropic energy (-Ex,iso; E^ iso ). The break time ranges from ~ hundreds 
to ~ 10 6 s with a typical value ~ 10 4 s, and the X-ray and the optical distributions are 
similar. The peak distribution of X-ray fluence is ~ 10 -6 ergcm -2 and in the optical band it 
is ~ 10 _8 ergcm~ 2 . The isotropic break luminosity (Lx,b or Er^) with a wider distribution, 
ranges from ~ 10 43 erg s _1 to ~ 10 49 erg s _1 for the X-rays with a typical value of 10 4 ' erg 
s^ 4 and from 10 42 erg s 1 to 10 4 ' erg s 1 for the optical band with a typical value of 
10 46 erg sW The isotropic energy release in the X-ray band (Ex,iso) generally ranges from 
~ 10 49 to ~ 10 53 erg, while for the optical band (£fu so ) it is from 10 48 to 10 54 erg. The peak 
distribution value of the X-ray band is ~ 10 52 erg and of the optical band ~ 10 5 °erg. 

Figure O shows the correlation between the X-ray and the optical properties for Group 
I. We examine the pair correlations for Aaq — Aa 2 and A a x — Aa 0 , and obtain r=0.40 with 
p=0.42 for Aaq — Act 2 and r=0.40 with p=0.55 for A a x — Ao 0 , which indicates that these 
are only weak correlations. Notices that the peak value of the distribution for T x ,b/E 0) b is 
around 0 in the logarithmic space. This indicates that a good fraction of Group I GRBs have 
the same break time in both the X-ray and the optical band. This is generally related to the 
achromatic break behavior. A detailed analysis of the chromaticity of afterglow behavior is 
presented in Wang et al. (2015). 

We compare the start (T Xjl , T o l ) and end times (T Xj2 , T o2 ) of the X-ray and optical band 
observations in Figure El The T x l — T o l panel (Figure El (a)) shows that the data points 
generally cluster around the equality line in the logarithmic space, and the X-rays range 
from ~ 10 2 seconds to ~ 10 5 seconds and the optical band from ~ 10 seconds to ~ 10 5 
seconds. Most of the data points are below the equality line in the logarithmic space for the 
T Xj 2 — T 0)2 panel (Figure El (b)), and the distribution is from ~ 10 4 to ~ 10' seconds for the 
X-ray band and from ~ 10 2 to ~ 10' seconds for the optical band. This indicates that the 
clear afterglow component was observed for the optical band earlier than for the X-ray band 
(it is likely that the early X-ray afterglow is usually polluted by the steep decay), and also 
for the end time of the observation. 


4. Comparison of the X-ray and Optical Temporal and Spectral Indices 

Our discussion concerns the slow cooling case, since most GRBs in our sample have the 
starting time (T Xj i) hundreds of seconds after the trigger. If the cooling frequency (u c ) lies 
between that of the X-ray (v x ) and the optical (u 0 ) bands (v m < u 0 < u c < u x ), then the 
differences in the temporal indices between the X-ray and the optical bands are the following: 


10 


(i) for the cases with no energy injection (q — 1), one has A a = | for the interstellar (ISM) 
medium scenario (Sari et ah 1998) and Aa = for the wind scenario (Chevaler & Li 
2000); (ii) for the energy injection case with q — 0 due to pulsar spin-down (Dai & Lu 1998; 
Zhang & Meszaros 2001), one has Aa = \ for the ISM case and Aa = — ~ for the wind 
case, which have the maximum difference e between the X-ray and optical bands (e.g. Oates 
et al. 2009); (iii) The difference in the spectral indices between the X-ray and the optical 
bands is at most A/3 = | for both the ISM and the wind models (the cooling break). If the 
cooling frequency is above the X-rays band ( v m < v Q < v x < is c ), one has the same spectral 
and temporal indices in both the X-ray and optical bands, regardless of the medium type 
and whether or not there is energy injection. The ranges of relations in the temporal and 
spectral indices are summarized in Table 5 (see also Table 2 of Zhang et al. 2006). The 
expressions of Aa and A/3 can be summarized in the following two formulae. 


Aa 


§, v m < v 0 < v c < u x , injection (q=0), (ISM) 

|, u m < u Q < v c < u x , no injection (q=l), (ISM) 

0, u m < v 0 < v x < u c , injection (q=0) and no injection (q=l), (ISM and wind) (10) 

I'm < Vo < V c < v x , no injection (q=l), (wind) 
v m < v Q < v c < v x , injection (q=0), (wind) 


A/3 


0, 

i 

2 ’ 


v m < v 0 < v x < v c , injection (q=0) and no injection (q=l), (ISM and wind) 
v m < v a < v c < v x , injection (q=0) and no injection (q=l), (ISM and wind) 


(ii) 


4.1. Temporal and spectral indices for the 4 groups 

In Figure [TQ1 we present the correlations between the temporal decay indices in the 
optical band versus those in the X-ray band (a x — a 0 panel). We present the correlation 
separately for our 4 groups. In the figure we mark the range deviating by -0.25 and 0.25 from 
the equality line, which is the range expected from the standard forward shock model for 
different circum-stellar media (ISM and Wind-like), as well as the range of -0.5, 0.5 (dashed 
lines) which is related to the energy injection model for different circum-stellar media (ISM 
and Wind-like). Note that if the points lying above the equality line decay more quickly in 
the optical than in the X-ray band, then the data are consistent with the wind-like model; if 
the points below the equality line decay more quickly in the X-ray than in the optical band, 
then the data are consistent with the ISM model. In the Figure fill we similarly present the 
correlation between the spectral indices in the optical band versus those in the X-ray band, 
for each of the 4 groups. We indicate the line corresponding to @ x — /3 0 lying within the range 
[0.0, 0.5] which is expected for different spectral regimes of the external shock model. 




-11 - 


The difference in temporal indices between the X-ray and the optical bands are compared 
to the empirical relations for various external shock models and are summarized in Table 4. 
For Group I (Figure [TUI (a)). 18 of 24 and 16 of 24 data points fall within the external shock 
model ranges (-0.5,0.5) for a X)1 — a o l and a Xi2 — a 0 , 2 , respectively. Out of these, 13 of 24 and 
10 of 24 data points fall within the standard afterglow model ranges (-0.25,0.25). Moving to 
Group IV (Figure fT0fd)h there are 27 of 30 data points which fall within (-0.5,0.5), which 
includes 13 of 30 data points within (-0.25,0.25) for a xj2 — a 0j2 . This indicates that the data 
for Group I and Group IV satisfy the external shock model for the global behavior. However, 
for Group II, although there are 17 of 24 data points which fall within [-0.5,0.5] and includes 

5 of 24 points within [-0.25,0.25] for ct Xil — a 0i2 , only 8 of 24 data points fall within (-0.5,0.5) 
and only rare points (only one case) within the ranges of (-0.25, 0,25) for cc Xi2 — a 0j2 . Similar 
to Group II, for Group III, 3 of 9 data points fall within [-0.5,0.5] and no point within [- 
0.25,0.25] for a Xi2 — a 0il . However there are 8 of 9 data points which fall within [-0.5,0.5] and 

6 of 9 points within [-0.25,0.25] for a Xi2 — a Q , 2 - These observations indicate that the data are 
generally inconsistent with the external shock model for the global behavior for Group II and 
Group III. Since these two groups include GRBs that potentially have a chromatic behavior 
(which is not predicted in the simplest afterglow model), this result is understandable: The 
GRBs that do not comply with the models may require two emission components to interpret 
the X-ray and optical emission, respectively. One of the two components should still come 
from the external shock. The data indeed show that for a certain section of the afterglow 
the data are still be consistent with the external shock model. A full investigation of this 
issue is presented in Wang et al. (2015). 

Based on the statistics of ct x — a 0 relation with the definitions of 4 Group samples, 
within the error, totally we find 54 of 87 (62%) GRBs that are consistent with the standard 
afterglow model (A a = [—0.25,0.25]). However, when more advanced modeling is invoked 
(e.g. long lasting central engine emission, reverse shock emission, or structured jets), 79 of 
87 (91%) GRBs may also satisfy the external shock models (Act = [—0.5, 0.5]). We also 
find that 53 of 87 (61%) GRBs correspond to the ISM model, and 34 of 87 (39%) GRBs 
correspond to the wind model. 

The spectral indices for the GRBs in each Group are shown in the /3 X — /3 0 panels of 
Figure [HI within error, the data generally fall well within the range of [0, 0.5], which is 
expected for the external shock models (the data of Group II have a little scatter around 
the external shock model ranges). Notice that we have collected the optical spectral indices 
from the literature. No time-resolved optical spectral analysis was available for most GRBs. 
A detailed time-resolved joint analysis of afterglow power density spectrum using X-ray and 
optical data will be presented in a separate paper (Li et al. 2015, in preparation). 
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In order to find the properties of the mean values of the temporal and spectral indices, 
in Figure [12] we show the distributions of the pre-, and post-break decay slopes, and the 
spectral indices for the X-ray and the optical bands and compare them0. All the distributions 
can be reasonably fitted with a Gaussian function in order to find the mean values. The decay 
slope of the pre-break slopes are 0.46±0.01 (optical) and 0.68±0.03 (X-rays). This indicates 
that a shallow decay segment is commonly detected in both the X-ray and optical band at 
the early time. The corresponding values of the post-break slopes are 1.23 ± 0.03 (optical) 
and 1.48 ± 0.02 (X-rays), for which the value of a Xi2 — a Qi 2 is 0.25, and fit very well with the 
case of an ISM medium, and slow cooling with spectrum regime (z/ m < v Q < v c < v x ). The 
spectral indices are 0.99 ± 0.23 (X-ray) and 0.69 ± 0.25 (optical), also well-consistent with 
the external shock model. 


4.2. Temporal indices for different epochs 

I 11 order to study evolution in the relations, we compare the X-ray and the optical 
temporal indices at different epochs, and divide the data into 4 epochs and test the data 
against the external shock model. 

We study a sample of GRBs that is selected according to the following: (i) We select 
the GRBs with redshift measurements (80 out of 87 bursts) so that their observed times can 
be corrected to the burst rest-frames; (ii) We chose four epochs in the rest frames: epoch 
1: < 10 3 s, epoch 2: 10 3 s — 10 4 s, epoch 3: 10 4 s — 10 5 s, and epoch 4: > 10 5 s. The reason to 
choose these different epochs is to investigate how the joint X-ray/optical emission properties 
are consistent or not consistent with each other in different epochs, (iii) We selected the 
cases for which there was at least one epoch during which both X-ray and optical data are 
available, (iv) If a light-curve break time (X-ray or optical) lies within one epoch, we include 
both the pre- and post-decay indices in this epoch. Using these selection criteria we finally 
get 50, 97, 80, 35 cases in the four epochs, respectively. For the different samples the number 
of cases are (19, 45, 41, 11) for Group I, (14, 26, 15, 7) for Group II, (5, 7, 5, 2) for Group 
III, and (12, 19, 19, 15) for Group IV. 

As shown in Figure IT3l in the a x — a Q panels for our epoch samples, consider the error, 
42 of 50 (84%), 76 of 97 (78%), 60 of 80 (75%) and 28 of 35 (80%) data points for Act fall 
within our external shock model ranges (-0.5,0.5) for epochs 1, 2, 3 and 4, respectively. Out 


7 The pre-break decay slope indices include 48 a x p (Group I+Group II) and 37 a 0 ,i (Group I+Group 
III). The post-break decay slope indices include 87 ct Xj 2 and 87 a 0 ,2 (for all Groups). The spectral indices 
include 87 /3x,2 and 80 do (for all Groups) 
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of these, 23 of 50 (46%), 45 of 97 (46%), 37 of 80 (46%) and 19 of 35 (54%) data points are 
within the standard afterglow model ranges (-0.25,0.25). This indicates that the data are 
generally consistent with the external shock models and there is no significant evolutionary 
effect in onr epoch samples. The results can be also seen in Table 3. 

In order to consolidate these results, and in particular to show that the data are not 
randomly distributed, we compare the observational data with a simulation using the Monte 
Carlo (MC) method. We simulated a random set of decay slopes with a Gaussian distribution 
in both the X-ray and optical bands, lying within the observed range of [0,3]. For each epoch 
the simulated set contains the same number of points as the observed set. The simulated 
points (red open circles) are also shown in Fig.13. We evaluate the consistency between the 
simulated and observed sets by a Kolmogorov-Smirnov (K-S) test and use the probability of 
the K-S test (Pk-s) to measure the consistency. A larger value of Pk-s indicates a better 
consistency. If a value of Pk-s is less than 0.1, then there is no consistency. Our results 
are shown in Table 5, which suggests that the simulated and observed data are generally 
inconsistent with each other. We conclude that the observed data sets are not randomly 
distributed. 


5. Confronting the Break Behavior with Various External Shock Models 

The canonical X-ray light-curve (Zhang et al 2006) can be described as an initially steep 
decay phase I (with a temporal power-law index, a x ~ —3 and F( t ) oc t a ), then changing 
to a shallow decay segment (or plateau) II (—1.0 < a x < 0). This phase is then followed 
by a normal decay segment III ( — 1.5 < a x < 1.0) and finally the light-curve changes into 
a post jet-break decay IV. Compared to the X-rays, the synthetic optical light-curves can 
be generally summarized as eight emission components with distinct physical original, the 
segments II-IV are usually observed as counterparts in optical afterglow light-curves (see 
Li et al. 2012, Liang et al. 2013). Thus, break behaviors in afterglow light-curves are 
commonly observed, the shallow-to-normal breaks at early times and the normal-to-jet break 
behavior at late times. Another possibility is a spectral break which is due to the cooling 
frequency crossing the X-ray or optical band. We summarized the properties of afterglow 
break behaviors, comparing the X-ray with optical bands and their statistic behaviors in 
Table 6. 


8 These eight components are: la. early prompt optical emission, lb. early reverse shock emission, II. 
early shallow decay segment, III. standard afterglow (normal decay), IV. post jet-break phase, V. optical 
flares, VI. late re-brightening emission, VII. late SN-GRB bump. 
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The external forward shock model provides well-predicted spectral and temporal prop¬ 
erties and allows for a wide variety of behaviors, it has been successful in explaining the 
broadband afterglow emission (see, e.g., Rees & Meszaros 1992; Meszaros & Rees 1993, 
1997; Paczynski & Rhoads 1993; Sari, Piran & Narayan 1998). The relations between spec¬ 
tral and temporal properties are denoted ’’closure relations”, which have been widely used 
to interpret the rich multi-wavelength afterglow observations. The temporal index a and 
spectral index [3 in various afterglow models and the typical numerical values of each case 
are summarized in Table 5 (see also Table 2 in Zhang et al. 2006; Sari et al. (1998, 1999); 
Chevalier & Li (1999)). 

The closure relations (a — f3) of the various external shock models in different regimes 
for all the cases are shown in Figure [£U together with the closure relations for the ISM and 
wind-like models (Liang et al. 2008). The X-ray data are shown in Figures [TT] (a) for the ISM 
model and [14] (b) for the wind-like model. For X-ray break samples (Group I+Group II), the 
pre-break segments (the black solid) are shallower than the model predictions according to 
the standard afterglow model, generally due to the energy injection causing a shallow decay 
segment. The post-break segments (the blue triangles) are well-consistent with the standard 
afterglow model, or fall well within the jet-break model region. We find the data for the X- 
ray non-break samples (the pink half-solid points, Group III+Group IV) are well-consistent 
with the standard afterglow model. The optical data are shown in Figure [L4] (c) and (d) 
for the ISM model and the wind-like modelj. It can be seen that optical data in jet-break 
regions (are clustered in our optical break samples: Group I and Group III) are fewer than 
X-ray data. By contrast, the fraction of GRBs which fall in the jet-break region for optical 
non-break samples is significantly larger than in the X-ray cases. A few GRBs fall in the 
jet-break regions with spectral regime I ( u opt > max{v m , u c )). 


5.1. Energy injection break 

The shallow decay segments were commonly detected in the early afterglow light-curves, 
generally believed to be caused by the GRBs outflow catching up with the external forward 
shock and with a long-lasting energy injection into the external shock which refreshes the 
afterglow emission. It can be defined by the afterglow theory, by which the temporal decay 
is shallower than what is predicted in the external shock model. 

We defined the shallow decay segment by the closure relation: an < 3 ^ 1 2 ~ 1 for the spectral 


9 Same conventions, but since there are no time-resolved optical spectral data analyses available, we use 
the same spectral index for both the pre- and post-break segments. 
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regime I (v > max^m^Vc)) for both the ISM and the wind-like models, «i < for ISM 
model and ag < ' 2 +1 for wind-like model for spectral regime II (u m < v < is c ), respectively. 
We find that generally a very high fraction of ag are in a shallow decay segment for both 
the X-ray and the optical break samples. As summarized in Table 6, 42 of 48 GRBs have a 
shallow decay segment for X-ray break samples (including: 21 GRBs for Group I and Group 
II, respectively). 27 of 30 GRBs with a shallow decay segment for optical break samples 
(including: 18 of 21 GRBs for Group I and 9 of 9 GRBs for Group III). No significant 
spectral evolution is observed for the shallow decay segment to the following phase for our 
42 X-ray, shallow decay samples (see also in Table 6), only 4 of 42 GRBs showed significant 
spectral evolution, indicating that the shallow decay should have an external-shock origin 
(see also Liang et al 2007). 

The energy injection behavior can be described as L(t) oc t~ q (see, e.g. , Zhang & 
Meszaros 2001), here q is the temporal index. The difference between the decay slopes 
before and after the time of the break depends on the spectral regime observed and the type 
of ambient medium, which can be summarized as in Table 5 (see also Zhang et al.(2006) in 
Table 2). This model is called the ’’energy injection model” (Zhang et al 2006) and produces 
a break in the light-curves, which are denoted ’’injection breaks” (Zhang et al. 2004 for 
review). The pre-break slope ag should give q = 0 for a constant luminosity, 


[ (q-l) + = (2p-6)+( P +% ; Um<u< i/ c (ISM) 

«i = { I + = - (2 P -2)+( P +% z/ m < Z/ < U C { wind) (12) 

[ (s- 2) + (2+^ = (2p—4)+(p+3)g ^ y > ^( ISM and wind), 

while the post-break slope ag should correspond to q =1 (for a constant energy fireball, the 
scaling law is the same as q— 1, Zhang & Meszaros 2001). Then, according to Zhang et al. 
(2006), one should have 


ag — 


3/3 __ 3Q~1) _ 3ai+3 
2 4 2 ’ 

3/3+1 _ 3p—1 _ 3 qi + 1 
2 4 2 

3/3—1 _ 3p—2 _ 3 qi+2 


I'm < V < ^c(ISM), 
z/ m < v < z/ c (wind), 
v > z/ c (ISM and wind), 


(13) 


The temporal and spectral properties of the afterglow after the break (the normal de¬ 
cay phase) should satisfy the ’closure relation’ of the external shock model (e.g. Zhang & 
Meszaros 2004), i.e. as described in Table 5. The relation of the a± — ag panel for all our 
break samples (Group I+II for the X-ray emission and Group I+III for the optical emission) 
is shown in Figure [T5l 
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Since the calculation of the q and p indices depend on which afterglow model is used, 
our first step is to use the observed quantity of the temporal index a and spectral index (3 to 
determine which spectral regime the burst belongs to. Spectral regime I ( v obs > max(u m , z/ c )) 
corresponds to the observed band lying above the cooling frequency; and Spectral regime 11 
(u m < z/ Q bs < r'c) corresponds the observed band lying below the cooling frequency. Since dif¬ 
ferent spectral regimes give different predicted decay indices through the “closure relations” 
of the external shock model (e.g., Table 2 of Zhang et al 2006; Gao et al. 2013), we can 
combine the a and (3 information to determine the spectral regime of a temporal segment 
in each GRB. Since the normal decay phase is not contaminated by energy injection (Zhang 
et al. 2007), we use its temporal index a to perform the analysis. We first use the temporal 
index a to calculate the theoretically expressed spectral index (ft) for the different spectral 
regimes: 

v > max(v m , z/ c )(ISM and wind) 

v m < v < z/ c (ISM) (14) 

v m < v < z/ c (wind) 

Second, we compare this index (j3' a ) to the observed spectral index (/3 0 b s ): X 2 — (P a ~ 
Aibs) 2 /( <7 ^(a)+ cr ^(obs))- We thereby can identify which spectral regime the observation belongs 
to. We find that for our X-ray observations, 19 of 87 GRBs belong to spectral regime I, and 
68 of 87 GRBs belong to spectral regime II. For optical emission the following was found: 
9 of 80 GRBs belong to spectral regime I, and 68 of 80 GRBs belong to spectral regime II. 
Note that the optical spectral index of 7 GRBs is not available and therefore we did not 
judge the spectral regime of the optical band for those GRBs. 

After determining the spectral regimes for each individual burst, the electron spectral 
index p can be determining through, 

v > max(iy m , z/ c )(ISM and wind) 

v m < v < z/ c (ISM) (15) 

u m < v < z/ c (wind) 

Similarly, the the luminosity injection index q (for energy injection cases, calculated 
by shallow decay segments) can be determined from the temporal index a and the spectral 
index /?, 


P = 



P a = ¥ 


2a+l 
3 ’ 
2 a 
3 ’ 
2a—1 


Q = 


201 - 2 / 31+2 

l+/3i ’ 

2ai-2/3i+2 

2+0i ’ 

2ai—2/3i 
1+01 ’ 


v > max(u m , v c ) (ISM and wind) 
v m < v < z/ c (ISM) 
v m < v < i/ c (wind) 


( 16 ) 
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The distributions of the p and q indices for both the X-ray and optical bands, are shown 
in Figure fl6l (a) and (b). We calculated the q index for onr break samples with the temporal 
index of pre-break in a shallow decay segment. The q value ranges within [-0.90, 0.88] for 
the X-ray band and within [-0.79, 0.90] for the optical band, and the best Gaussian function 
fitting with a center value g x ( c ) = 0.46 ± 0.03 and g Q ( c ) = 0.13 ± 0.11 for the X-ray and the 
optical band, respectively. The p index is derived from the observed spectral index of the 
post-break decay slope for all the cases. The p value ranges within [1.78, 4.18] and 6 of 87 
cases had p < 2.0 for the X-ray band and within [1.28, 4.42] and 17 of 87 cases has p < 2.0 
for the optical band. The best Gaussian fitting with a center value p x ( c ) = 2.63 ± 0.02 and 
Po(c) — 2.58 ± 0.04 for the X-ray and the optical band, respectively. 

The A ct for each case is summarized in Table 4. The functional dependence of A a—p, q 
for the different external shock model cases can be summarized as: 

( \{p + 2)(1 - q ), v > z/ c (ISM and wind), 

Ace = < \{p + 3)(l-q), u m <u< z/ c (ISM), (17) 

l I (P + 1)(1 - q), Vm< v < v c (wind), 

The two-dimensional panels between the observed quantity Ace and the value of the model 
parameters p and q are shown in Figure [T6] (c). For the Ace — p panel, the various external 
shock model lines are well covered by all the data points when q=- 0.5, 0, 0.5, respectively. 
For the Ace — q panel, we also find the external shock model lines are well covered by the 
data, when we have p= 2.0, 3.0, respectively. Compared with the optical band, we find that 
the X-ray data are gathered in the regions which have higher p values. 


5.2. Jet break 

Another common break behavior in the afterglow is a post-jet break. Physically there 
are two phenomena that could cause a steepening decay in the late afterglow light-curves: 
the edge effect and sideways expansion. Act = 0.75 (ISM) and Act = 0.5 (wind-like) are 
expected from the edge effect. A sharper break (from a norma i ~ —1.2 to aj et ~ —p) is 
suggested by sideways expansion, the post-jet break decay index is approximately ~ —p. 
Thus the criterion for a good candidate for jet-break can generally be summarized as: (i) 
ct 2 ~ —p. Therefore, we take ct 2 > 1.5 as our definition of a jet-break sample (see also 
Liang et al. 2008), since this value is difficult to interpret by the standard fireball model, 
(ii) Act ranges within [0.5, 1.2] for normal-to-jet breahl 10 l. (iii) The ’Closure relation’ should 


10 Here the Act=1.2 is supposed to be a typical value of p. 
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be satisfied. As shown in the Figure [141 in the a — (3 panel, a good fraction of GRBs in 
our samples (both the X-ray and optical band) are consistent with the expectations of the 
external shock jet models, especially in the wind-like model, (iv) No spectral evolution. 
Since jet-break behavior is purely of dynamical effects, there should be no spectral evolution 
across the break and achromatic behavior should be observed in the multiple-band afterglow 
light-curves. 

The distribution of the p index is shown in Figure dHl and the peak distribution of the p 
value is generally 2.5 for both the X-ray and the optical bands. As summarized in our Table 
6, there are no jet-break cases with a spectral evolution. 32 of 48 GRBs for the X-ray break 
cases (including: 14 of 24 for Group I and 18 of 24 GRBs for Group II, respectively) and 13 
of 33 GRBs for optical break cases (including: 11 of 24 GRBs for Group I and 2 of 9 GRBs 
for Group III) satisfy the criterion of > 1.5. Based on this criterion, the A a x ranges 
within [0.56, 2.27] for the X-ray and A a 0 within [0.75, 1.84] for the optical jet-break cases, 
respectively. This indicates that the data are difficult to fully interpret by the criterion for 
a normal-to-jet break decay segment (Aa ranges within [0.5, 1.2] expected by the theory), 
since we find that there are also 16 of 32 GRBs for the X-ray and 10 of 13 GRBs for the 
optical band, with A a > 1.2. We find that the typical time for our jet-break samples is 
~ 10 5 seconds, and for our energy injection break samples ~ 10 1 seconds. However, it is 
interesting that we find that aq is generally in a shallow decay segment and aq is generally 
in the post-jet break segment. All these cases have Aa > 1.2, indicating that the energy 
injection break time appears similar to the jet-break times for those GRBs. This is the 
reason for Act being large. 


5.3. Spectral break 

Another possibility that can cause a break in the X-ray (or optical) light-curve is due 
to spectral evolution where the spectral peak of the synchrotron emission moves through the 
X-ray band. However, since, in general, there is a lack of significant spectral evolution, this 
possibility is largely excluded (De Pasquale et ah 2006; Nousek et ah 2006). If the cooling 
frequency ( v c ) crosses the X-ray band, the change of the X-ray spectrum should be reflected 
by A f3 x = | for the ISM scenario and A /3 X = — | for the wind-like circum-burst medium 
scenario (Sari et ah (1999); Chevalier & Li (1999)). The change in the light-curve should 
then correspond to A a x = -)■ for both the ISM and wind-like scenario. 

We show the distributions of Aa for all our cases with breaks in Figure [171 We have 44 
X-ray break cases (Group I+II), and the distribution of A a x is in the range [0.44,2.27], which 
is larger than A a = 0.25 as expected by the cooling frequency crossing the X-ray band model. 
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This is similar for the 33 cases with optical afterglow breaks (Group I+11I). The distribution 
of A a Q is in the range of [0.47, 1.84]. All bursts have A a values that are larger than the 
value expected by the cooling frequency crossing the optical band (Aq = 0.25). 

On the other hand, we derived the photon spectral indices of the pre-break and post¬ 
break portions of the X-ray afterglow for our break samples (Group 1+ If) and then calculated 
the spectral indices using /3 X = T x — 1. The distribution of the spectral indices of the pre¬ 
break and post-break portions are shown in Figure [171 (b). and the center values with fitting 
by Gaussian functions to the distributions are /? Xi i = 0.95 ± 0.03 and = 0.98 ± 0.01. 
The /3 X) 2 — /3 x ,i panel shows that all the bursts (to within the error) are in the vicinity of the 
equality line. This indicates that the pre-break and post-break portions of the light-curves 
lack the significant spectral evolution in most of the cases. Since no optical time-resolution 
was available for analysis of all the optical break sample, we are not able to examine the 
change of the spectral index for the transition from pre-break to post-break. 

Both the X-ray and optical afterglow break behavior can thus not be interpreted by 
the cooling frequency crossing the X-ray or optical band. Indeed, there is also a lack of any 
significant spectral evolution from pre-break to post-break. 


6. Conclusions and Discussion 

We extracted the cleanest possible signals of the afterglow emission components, and 
include the observations not only by Swift/UVOT, but also by other ground-based telescopes, 
and obtain a sample of, in total, 87 bursts that have good quality, simultaneous detections in 
both the X-ray and the optical bands. We divided the sample into four groups depending on 
the existence of breaks in their afterglow light-curves, and examine their behavior in light of 
various external shock model empirical relations. We find that the data generally satisfy the 
various external shock empirical relations. Our main results can be summarized as follows: 

• About 10% Swift-GRBs (until the end of 2013) had good simultaneous observational 
data in both the X-ray and optical afterglow. Studying this sample in detail, we find 
that the X-ray and optical afterglow light-curves show similar statistical properties. 

• Characterizing the light-curves in 4 different groups (Fig.l), we investigate how well 
GRBs in each group satisfy the external shock model predictions. We find that in 
groups I and IV, a majority of GRBs follow the predictions the external shock model. 
The bursts in groups II and III have, per definition, chromatic light-curves, so that one 
cannot interpret the full evolution directly with the external shock model. However, a 
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large fraction of these burst have individual segments in the light-curves that do satisfy 
the relations of the external shock model. This is consistent with the following picture: 
Four these bursts, one cannot interpret both X-ray and optical data simultaneously 
in the external shock model. Only one band may be of the external shock origin. 
Part of the other band may arises from a different emission component (e.g. long 
lasting central engine emission, reverse shock emission, or structured jets) causing the 
chromatic behaviour, see further discussion in Wang et al. (2015). 

• Overall we find 62% GRBs that are consistent with the standard afterglow model 
(Act = [—0.25, 0.25]), when more advanced modeling is invoked, up to 91% of the bursts 
in our sample may be consistent with the external shock model (Act = [—0.5, 0.5]). 
Of these, 61% correspond to a ISM medium and 39% GRBs correspond to a wind 
medium. 

• The fraction of the bursts that are consistent with the external shock model is inde¬ 
pendent of the observational epochs (measured in the rest frame of GRBs). 

• Approximately half of all bursts have light-curve for which the break time of the 
shallow-to-normal decay is similar to that of the normal-to-jet decay, making an imme¬ 
diate change from the plateau phase to the jet phase in both the X-ray and the optical. 
We suggest that the normal decay segment for these bursts is short and therefore not 
detected. 

• The average values of the post-break decay-slopes, in the X-ray and optical band, are 
consistent with the case of an ISM medium and slow cooling electrons (u m < u 0 < z/ x < 

^c) • 

• We do not detect one single case in our sample with a temporal break that can be 
interpreted as the cooling frequency crossing the X-ray or optical bands. 

The models we have tested are the simplest analytical external shock models of GRBs. 
More sophisticated afterglow modeling has been carried out, which gives somewhat different 
predictions than the simplest models. For example, the recent numerical simulations of 
Duffcll & MacFadyen (2014) simulated hydrodynamical evolution of a jet starting from a 
collapsar engine to the late afterglow phase. The change of decay slope in their simulated 
light curve from the plateau to the following power-law decay is A a x = 9/8. We investigated 
48 GRBs in X-rays and 33 GRBs in optical that show such a break, and find that 23 of 48 
(48%) GRBs for the X-ray and 17 of 33 (52%) for the optical are consistent with having 
A a ~ 9/8, which is consistent with their simulation results. More detailed analyses are 
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needed to fully confront this and other more sophisticated modeling with the observational 
data. 

One important question is whether the X-ray and optical emissions come from the 
same emission component, as predicted by the standard afterglow model. One interesting 
fact is that we have identified 9 of 87 GRBs, whose afterglow light-curves show a genuine 
afterglow behavior, namely a single power-law decay sometimes with an early afterglow onset 
behavior, for both the X-ray and the optical emission. These GRBs are fully consistent with 
the simplest standard fireball model. There are no flares, no shallow decay segments, no 
re-brightening at late time (without extra energy injection), no post-jet breaks. The nine 
candidates are: GRBs 050603, 061007, 080804, 090323, 090328, 090902B, 090926A, 120815A, 
130427A. 

Other GRBs have more complicated features. Testing whether these GRBs are con¬ 
sistent with the external shock models requires more detailed modeling. In this paper, we 
performed a series of consistency checks, including the a — [3 closure relations in individual 
segments and epochs and Act and A/3 predictions. Onr results show that in general most 
of them are still consistent with the external shock model predictions. Nonetheless, in order 
to address whether these GRBs are fully consistent with the afterglow models, much more 
detailed analysis and modeling, including the achromaticity and global consistency with the 
closure relations, are needed. This has been carried out in a separate work (Wang et al. 
2015). That work also suggests that most GRBs are consistent with the external shock 
model of GRBs, which is fully consistent with the finding of this paper. 
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Table 1. XRT Observations and Fitting Results 
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050502A 

2.87,1050 

18.8=1=2.5 

0.66=1=0.06 

1.66=1=0.05 



0=1=29 



050824 

6.31,2030 

74.9=1=58.4 

0.32=1=0.14 

0.94=1=0.13 

30=1=17 

2.19=1=1.0 

5±3 

1.85=1=0.19 

2.48=1=0.38 

051221A 

6.03,848 

26.7=1=5.2 

0.18=1=0.15 

1.43=1=0.06 

12=1=3 

1.70=1=0.32 

1±0 

1.87=1=0.22 

1.97=1=0.16 

060111B 

0.20,222 

7.6=1=4.0 

0.86=1=0.07 

1.58=1=0.14 

20=1=9 



2.03=1=0.2 

2.07=1=0.32 

060418 

0.36,537 

1.3=1=0.4 

0.87=1=0.13 

1.57=1=0.05 

47=1=19 

1390=1=590 

26=1=11 

1.98=1=0.22 

1.95=1=0.22 

060906 

0.20,258 

12.6=1=1884.1 

0.29=1=0.06 

1.82=1=0.17 

12=1=2 

343=1=51 

33=1=6 

2.07=1=0.21 

1.77=1=0.25 

060908 

0.08,450 

0.7=1=0.2 

0.44=1=0.10 

1.56=1=0.07 

36±9 

3360=1=920 

31±8 

2.05=1=0.29 

1.98=1=0.22 

060927 

0.09,103 

4.2=1=2212.5 

0.68=1=0.12 

1.99=1=0.46 

7=1=5 

1340=1=980 

36=1=24 

1.86=1=0.25 

1.96:1=0.35 

061121 

0.30,2100 

6.7=1=0.5 

0.41=1=0.03 

1.48=1=0.02 

261=1=27 

849=1=79 

116=1=12 

1.95=1=0.12 

1.82=1=0.06 

070318 

5.59,708 

341.2=1=62.2 

0.85=1=0.03 

2.41=1=0.51 

40=1=8 

0.818=1=0.17 

7=1=1 

2.02=1=0.12 

1.82=1=0.49 

070420 

0.25,525 

3.4=1=0.3 

0.26=1=0.06 

1.47=1=0.02 

113=1=11 



2.07=1=0.15 

1.92=1=0.12 

071031 

6.13,344 

104.2 

0.97=1=0.08 

1.75=1=0.48 

4=1=1 

4.68=1=0.69 

6=1=1 

1.73=1=0.3 

1.26=1=0.98 

080810 

3.85,366 

17.2=1=9.6 

1.17=1=0.14 

1.87=1=0.11 

32=1=22 

525=1=440 

73±52 

2.06=1=0.12 

2.2=1=0.22 

081008 

0.52,261 

18.5=1=4.7 

0.86=1=0.04 

1.85=1=0.14 

60=1=20 

169=1=52 

56=1=19 

2.05=1=0.13 

1.83±0.32 

081126 

0.07,368 

5.8=1=1.2 

0.60=1=0.04 

1.56=1=0.06 

85=1=15 



1.89=1=0.17 

2.28=1=0.19 

081203A 

0.14,345 

14.9=1=2.8 

1.16=1=0.01 

2.14=1=0.12 

87=1=27 

190±50 

92=1=28 

2.02=1=0.11 

2.08=1=0.35 

090510 

0.10,67 

1.4=1=0.2 

0.62=1=0.04 

2.15=1=0.09 

57=1=9 

459=1=70 

12=1=2 

1.67=1=0.1 

2.05=t0.2 

100418A 

0.40,2090 

125.0=1=17089.6 

-0.OliO.05 

1.66=1=0.10 

49=1=8 

2.31=1=0.31 

5=1=1 

1.89=1=0.29 

1.84=1=0.17 

100621A 

0.40,1770 

99.3=1=20.6 

0.86=1=0.03 

1.66=1=0.08 

371=1=51 

6.37=1=1.5 

28=1=4 

2.34=1=0.1 

3.06=1=0.31 

100814A 

3.88,4850 

148.9=1=7.33 

0.50=1=0.02 

2.09=1=0.06 

215=1=15 

51.52.6 

113=1=8 

1.896=1=0.048 

2.08=1=0.13 

100906A 

0.30,202 

13.2=1=1.3 

0.76=1=0.03 

2.12=1=0.80 

172=1=28 

560=1=80 

127=1=20 

1.96=1=0.1 

2.05=1=0.15 

120119A 

0.20,317 

32.5=1=7.0 

1.00=1=0.02 

2.03=1=0.18 

86=1=26 

66.7=1=17 

64=1=19 

1.61=1=0.11 

1.92=1=0.38 

121217A 

5.55,1340 

24.9=1=4.5 

0.47=1=0.09 

1.43=1=0.05 

134=1=22 



2.06=1=0.16 

2.05=1=0.16 

Group III 

050416A 

1.00,5800 



0.88=1=0.01 

47=1=6 


5=1=1 


2.27=1=0.1 

050525A 

5.92,1080 



1.54=1=0.04 

32=1=10 


3=1=1 


2.09=1=0 

050801 

0.19,442 



1.07=1=0.06 

6=h3 


4=1=2 


1.86=1=0.21 

051109A 

3.56,1520 



1.20=1=0.01 

96=1=11 


122=1=15 


2.05=1=0.07 

051111 

5.63,48 



1.61=1=0.06 

6±3 


4±2 


2.23=1=0.17 

061126 

1.61,2460 



1.34=1=0.01 

151=1=14 


53=1=5 


1.93=1=0.07 

080413A 

0.14,13 



1.21=1=0.03 

8=1=2 


11=1=3 


2.15=1=0.24 

090426 

0.13,497 



1.02=1=0.04 

8=1=2 


12=1=3 


2.12=1=0.17 

130702A 

89.14,15100 



1.16=1=0.01 

418=1=71 


2=1=0 


1.94=1=0.08 

Group IV 

050603 

34.81,1050 



1.71=1=0.06 

15=1=9 


26=1=15 


1.97=1=0.13 

050721 

0.20,136 



1.32=1=0.02 

71=1=5 




1.89±0.19 

050820A 

4.64,3680 



1.24=1=0.01 

236=1=28 


361±44 


1.98=1=0.05 

060124 

11.44,1960 



1.33=1=0.01 

202=1=30 


248=1=37 


1.98=1=0.05 

060206 

10.00,3430 



1.25=1=0.02 

50=1=10 


156=1=32 


1.9=1=0.19 

061007 

0.09,407 



1.73=1=0.01 

471=1=21 


193±9 


2.02=1=0.09 






















Table 1—Continued 


GRB 

Txl,T x 2 a 

T x , b a 

a= x ,i a c* x ,2 a 

Sx b 

L X,b° 

E X,lso C 

r c 

1 x i 

r c 

1 x 2 

070125 

46.83,1020 


1.63±0.12 

32=1=46 


19±28 


2.06±0.15 

070311 

7.10,54 


1.09±0.06 

10±4 





071010A 

34.51,472 


1.89±0.14 

20±2 


5±1 


1.85±0.21 

071025 

0.15,257 


1.54±0.01 

113±10 


462±39 


2.24±0.13 

071112C 

0.09,170 


1.33±0.01 

41±5 


7±1 


1.7=t0.13 

080319B 

0.06,2730 


1.61±0.00 

3001±68 


689±16 


1.91±0.07 

080804 

0.11,426 


1 . 11 ± 0.01 

31±3 


35±3 


1.82±0.1 

080928 

4.28,222 


1.59±0.03 

20±6 


14±4 


2.14±0.1 

081109A 

0.33,477 


1.22±0.02 

47±6 


12±1 


2.2=t0.13 

090102 

0.40,696 


1.30±0.01 

110±8 


66±5 


1.8=t0.08 

090323 

70.13,1050 


1.62±0.09 

11±9 


28=t23 


2.2=t0.28 

090328 

57.28,909 


1.68±0.08 

26±17 


4±2 


2.09±0.26 

090902B 

45.21,1460 


1.40±0.03 

49±17 


40=t 14 


1.82±0.12 

090926A 

46.68,1810 


1.41±0.04 

32±15 


34±16 


2.15±0.13 

091029 

0.70,1830 


1.27±0.05 

39±9 


65±16 


2.16±0.09 

100728B 

3.55,105 


1.51±0.07 

6±2 


6±2 


2.08±0.18 

100901A 

8.00,1560 


1.48±0.03 

77±6 


39±3 


2.29±0.09 

101024A 

4.28,137 


1.37±0.04 

12±5 




1.72±0.16 

110205A 

0.73,246 


1.67±0.01 

54±6 


62±7 


1.95±0.08 

110213A 

0.15,485 


2.04±0.04 

165±18 


89=t 10 


2.00±0.06 

110918A 

107.45,3880 


1.61±0.03 

302±73 


76=t 19 


1.96±0.09 

120404A 

0.80,21 


1.79±0.12 

16±4 


28±7 


1.91±0.13 

120815A 

2.73,41 


0.92±0.04 

19±7 


24±9 


1.83±0.12 

130427A 

0.35,15700 


1.28=h0.00 

4837±98 


139±3 


1.79±0.04 


a The time intervals, for the fit results of the X-ray afterglow for each GRB; and the break time for the X-ray 
break cases, in units of ks; the decay slopes before and after break time. 

^X-ray fluence, which is calculated by integrating the flux over the time intervals, in units of 10 — ® erg cm - 
Note that we corrected the X-ray fluence from absorbed to unabsorbed by calculating the ratios of the absorbed 
and unabsorbed flux, which collect from the Swift archive. 

c The break luminosity in the rest frame, in units of 10^ erg s~ ^; and the X-ray isotropic total energy release, in 
units of 10 50 erg; the photon spectral index of X-ray before and after break time. 



















Table 2. Optical Afterglow Observations and Fitting Results 


GRB(Band) 

T 0 i, T Q 2 a 

T o,b a 

<*o,l b 

«o,2 b 

Sr c 

L R,b° 

SR,iso C 

0o d 

2 d 


Group I 

050319(R) 

0.04,994 

249.6il3.72 

0.54i0.004 

1.94 

304i21 

2.81i0.13 

667i46 

0.74i0.42 

3.24^ 

050408(R) 

3.35,215 

33.8i2.33 

0.46 

1.12 

30i2 

0.22i0.01 

12il 

0.28i0.33 

1.2357 <2) 

050730(R) 

0.07,73 

8.5i0.66 

0.47i0.02 

1.45i0.03 

276i28 

117.34i7.49 

839i84 

0.52i0.05 

3.969 (1) 

050922C(R) 

0.74,606 

6.4=1=0.85 

0.46i0.03 

1.50i0.02 

386i73 

46.41i6.19 

439i83 

0.51i0.05 

2.198 <3) 

060210(R) 

0.06,7 

1.1=1=1.08 

0.11i0.40 

1.22i0.13 

18i23 

82.24i98.43 

55i70 

0.37i0.1 

3.91 (3) 

060526(R) 

0.06,894 

124.2=1=2.74 

0.73i0.01 

2.24i0.02 

535i24 

6.81i0.20 

1156i52 

0.51i0.32 

3.2l( 3 > 

060605(R) 

1.76,112 

25.0=1=2.03 

0.94i0.04 

2.59i0.13 

240i49 

32.66i5.44 

675il39 

1.06 

3.78a) 

060614(R) 

6.16,1280 

66.3=1=0.91 

0.16i0.03 

2.OliO.02 

283i9 

(100.00±1.00)xl0~ 4 

liO 

0.47i0.04 

0.125 <3) 

060708(R) 

0.06,208 

0.7=1=0.21 

0.06i0.10 

0.85i0.03 

358il05 

96.09il8.37 

320i94 


0.249 <4) 

060714(R) 

3.86,185 

7.7 

0.20 

1.17i0.05 

93i6 

20.20i0.40 

151i9 

0.44i0.04 

2.7l( 3 ) 

060729(U) 

0.70,662 

53.4=1=6.12 

0.16i0.10 

1.39i0.06 

3946i646 

2.43i0.24 

296i48 

0.78i0.03 

0.54a) 

060912A(R) 

0.10,24 

0.2(fixed) 

0.78 

1.32i0.02 

62i2 

46.61il.15 

14i0 

0.62 

0.8 (1 ) 

061021(R) 

0.09,298 

88.8=1=13.25 

0.65i0.01 

2.OliO.24 

330i66 

0.03i0.00 

10i2 


0.3463 <5 ) 

070411(R) 

0.70,517 

286.6=1=4.17 

0.95i0.00 

2.27 

197i6 

0.70i0.02 

371ill 


2.954a) 

070518(R) 

2.11,312 

30.0±38.86 

0.60i0.12 

1.80i0.12 

14i21 

O.lOiO.Ol 

5i7 

0.8 

1.16 (5 ) 

080310(R) 

0.30,124 

2.1=1=0.19 

0i0.03 

1.12i0.02 

198il8 

88.59i4.08 

267i24 

0.42i0.12 

2.43a) 

080413B(R) 

0.08,1901 

372.7=1=7.77 

0.50i0.01 

3.07 

400i50 

0.19i0.01 

126il6 

0.25i0.07 

l.g 7 ) 

080710(R) 

2.54,267 

8.7=1=0.08 

0.45i0.004 

1.52i0.01 

382i7 

4.93i0.05 

71il 

0.8i0.09 

0.85a) 

080721(R) 

0.16,2640 

20.0(fixed) 

1.OliO.01 

1.50 

460i9 

13.41i0.30 

700il3 

0.68i0.02 

2.602 <8 ) 

081029(R) 

4.95,253 

16.9=1=2.80 

0.48i0.19 

2.08i0.09 

326i98 

125.76i26.15 

943i283 

1.08i0.02 

3.847^) 

090618(R) 

0.08,73 

30.5=1=1.46 

0.70i0.002 

1.51i0.06 

836i51 

0.70i0.03 

63i4 

0.5 

0.54ao) 

090727(R) 

0.47,52 

3.3=1=3.10 

0.07i0.33 

0.78i0.16 

15i9 

0.90i0.41 

7i4 

0.81i0.28 

0.971a 1 ) 

091018(R) 

0.15,273 

64.1=1=1.79 

0.90i0.00 

1.67 

468i20 

0.49i0.01 

115i5 

0.58i0.07 

091127(1) 

7.93,531 

27.5=1=7.73 

0.45i0.09 

1.48i0.13 

774i310 

0.86i0.20 

48il9 

0.18 

0.49a 2 ) 

Group II 

050401(R) 

0.04,1120 



0.75i0.04 

33il6 


60i29 

0.39i0.05 

2.9 W 

050502A(V) 

0.05,18 



l.lli0.02 

115il7 


326i49 

0.76i0.16 

3.793a) 

050824(R) 

0.63,8990 



0.56 

725il3 


131i2 

0.4i0.04 

0.83a) 

051221A(R) 

11.12,445 



0.95i0.03 

30ill 


2il 

0.64i0.05 

0.5465a) 

060111B(R) 

0.07,13700 



1.12i0.08 

30i20 



0.7 

1.489a) 

060418(H) 

0.08,8 



1.27i0.02 

469i28 


263il6 

0.78i0.09 

060906(R) 

0.66,5 



0.83i0.17 

9i2 


23i4 

0.56i0.02 

3.685a) 

060908(R) 

0.06,1110 



1.09i0.01 

112il0 


97i9 

0.3i0.06 

1.8836a) 

060927(V) 

0.77,1 



1.17i0.04 

9i4 


49i20 

0.860.03 

5.6a) 

061121(R) 

0.05,1 



1.07i0.23 

165il7 


73i7 

0.95 

1.314a 2 ) 

070318(V) 

0.06,87 



1.14i0.06 

263i42 


48i8 

0.78 

0.836 (7 ) 

070420(R) 

0.12,11 



0.90i0.08 

167i33 




2.69a) 

071031(R) 

0.29,26 



0.84i0.001 

48i0 


77i0 

0.64i0.01 

080810(R) 

0.04,498 



1.20i0.01 

898il02 


2079i235 

0.44 

3.35a o ) 

081008(R) 

0.11,185 



0.96i0.00 

98i2 


92i2 

1.08i0.02 

1.967a 4 ) 

081126(R) 

0.10,1 



0.39i0.01 

49i0 




2.l(°) 

081203A(U) 

0.08,6 



1.47i0.00 

1244i6 


1306i7 

0.596 

090510(R) 

0.11,104 



0.98i0.12 

4i2 


liO 

0.76i0.14 

0.903^°) 

100418A(R) 

1.10,1370 



1.37i0.13 

158i80 


16i8 

0.98i0.09 

0.6235a 5 ) 

100621A(R) 

0.26,4 



0.46i0.06 

5il 



0.8i0.1 

0.542a°) 

100814A(R) 

0.52,18 



1.90i0.05 

28i35 



0.180.08 

1.44a 7 ) 

100906A(R) 

0.05,11 



1.07i0.02 

487i35 


246il7 


1.727a 8 ) 

120119A(R) 

0.08,17 



1.36i0.03 

69i6 


51i5 

0.89i0.01 

1.728a 9 ) 

121217A(R) 

1.34,347 



0.78i0.01 

147i3 



0.87i0.04 

3 1 ( 20 ) 

Group III 

050416A(R) 

4.09,33 

15.3=1=6.84 

0.39i0.15 

1.32i0.32 

7i4 

0.05i0.02 

liO 

1.3 

0.65 w 

050525A(V) 

0.07,316 

10.4=1=2.33 

0.78i0.01 

1.52i0.07 

202i61 

0.49i0.10 

19i6 

0.97i0.1 

0.606a) 

050801(R) 

0.02,22 

0.2=1=0.01 

0.04i0.02 

1.20i0.01 

131i5 

196.77i4.77 

80i3 

li0.16 

1.56a) 

051109A(R) 

0.04,1040 

0.4=1=0.32 

0.40i0.08 

1.04i0.04 

313il93 

210.95il09.42 

398i246 

0.7 

2.346a) 

051111(R) 

0.03,8 

3.1=1=0.67 

0.81i0.01 

2.10i0.70 

299i86 

31.38i5.75 

181i52 

0.76i0.07 

1.55a) 

061126(R) 

1.06,156 

30.4i8.94 

0.72i0.05 

1.39i0.07 

114i46 

0.63i0.20 

40il6 

0.95 

1.1588a) 

080413A(R) 

0.01,1 

0.4=1=0.07 

0.64i0.03 

1.82i0.29 

136i33 

477.21i88.04 

185i44 

0.67 

2.433a) 

090426(R) 

0.09,2 

0.2=1=0.05 

0.32i0.14 

1.24i0.06 

14i3 

113.06il7.79 

21i5 

0.76i0.14 

2.609a) 

130702A(R) 

15.13,280 

91.3i81.29 

0.44i0.22 

1.22i0.29 

1259il067 

0.03i0.02 

6i5 

0.7i0.1 

0.145a 1 ) 

Group IV 

050603(V) 

34.09,220 



1.80 

318i24 


556i42 

0.2i0.1 

2.821^ 

050721(R) 

1.48,8 



1.32i0.03 

26i7 



1.16i0.35 

2.612a) 

050820A(R) 

0.23,663 



0.91i0.00 

981il4 


1503i21 

0.72i0.03 

060124(R) 

3.34,1980 



0.88i0.02 

1015i247 


1246i303 

0.73i0.08 

2.296a) 

060206(R) 

18.41,202 



1.32i0.01 

339i43 


1061il35 

0.73i0.05 

4.408a) 

061007(R) 

0.03,15 



1.67 

2657i53 


1090i22 

0.78i0.02 

1.26a) 
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Table 2—Continued 


GRB(Band) 

Tol>^o2 a 

To,b a «o,l b a o,2 b 

Sr c 

^R,b° ^R.iso 0 

Po d 

2 d 

070125(R) 

105.86,349 

1.64=1=0.15 

349=1=1201 

210=1=724 

0.55=1=0.04 

0.82( 22 ) 

070311(R) 

0.07,92 

0.73=1=0.02 

341=1=86 


1=1=0.2 

0.98 (6) 

071010A(R) 

12.24,523 

1.91±0.02 

256±22 

64=1=5 

0.68 

0.98 (6) 

071025(J) 

0.17,15 

1.03=1=0.04 

131±11 

532=1=46 

0.42=1=0.08 

4.2< 6 ) 

071112C(R) 

0.13,70 

0.88=1=0.02 

52±6 

9=1=1 

0.63±0.29 

0.82( 23 > 

080319B(R) 

0.18,4590 

1.31±0.00 

1564=1=50 

359=1=11 

0.66 

0.94(®) 

080804(R) 

1.16,26 

0.93=1=0.02 

29=1=7 

33=1=8 

0.43 

2.2^ 6 ^ 

080928(R) 

7.77,233 

2.13=1=0.08 

187=1=138 

133=1=98 

1.08=1=0.02 

1.69 (1 0) 

081109A(H) 

0.17,67 

0.94=1=0.03 

86=1=18 

21=1=5 



090102(R) 

0.94,265 

1.02=1=0.01 

49=1=6 

30=1=4 

0.74 

1.547 (6) 

090323(R) 

162.43,769 

1.77=1=0.19 

75=1=1603 

192=1=4114 

0.65=1=0.13 

3.57 (1 °) 

090328(U) 

57.89,263 

1.54=1=0.12 

228=1=475 

32=1=67 

1.19=1=0.2 

0.736U°) 

090902B(R) 

62.29,564 

1.00=1=0.06 

98=1=74 

80=1=60 

0.68=1=0.11 

1.822 (24) 

090926A(R) 

84.92,102 

1.47±0.01 

74=1=8 

78=1=9 

0.72=1=0.17 

2.1062< 25 > 

091029(1) 

0.31,5 

0.59=1=0.01 

14=1=0 

24=1=0 

0.57=1=0.12 

2.752 (26) 

100728B(R) 

0.16,6 

1.01±0.02 

13=1=2 

14=1=2 


2.106 (27) 

100901A(R) 

10.19,543 

1.38=1=0.001 

949=1=30 

480=1=15 

0.82 

1.408 (28) 

101024A(R) 

1.42,160 

1.02=1=0.08 

22=1=15 




110205A(R) 

0.54,384 

1.47±0.00 

355=1=3 

410=1=3 

0.33 

2 22( 29 ) 

110213A(R) 

0.10,183 

0.87±0.08 

1515=1=421 

820=1=228 

1.12=1=0.24 

1.46( 29 ) 

110918A(R) 

122.43,805 

1.31±0.04 

434=1=252 

110=1=64 

0.7=1=0.02 

0.984< 30) 

120404A(R) 

0.73,57 

1.47±0.02 

177=1=6 

228=1=7 

1.05=1=0.09 

2.87( 31 ) 

120815A(R) 

0.17,11 

0.64±0.001 

40=1=1 

51=1=1 

0.78=1=0.01 

2.358< 32 > 

130427A(R) 

0.31,456 

1.04±0.001 

5855=1=33 

168=1=1 

0.69=1=0.01 

0.34^ 33 ^ 


References. — (1) Butler et al.(2007); (2) Nysewander et al.(2009); (3) Nava et al.(2008); (4) Jakobsson et al(2006GCN,5319) 
(5) Fynbo et al. (2009); (6) Butler et al.(2010) (7) Krimm et al.(2009) (8) Ukwatta et al. (2010); (9) Nardini et al. (2011) (10) 
Guetta et al.(2011); (11) Wiersema et al.(2012) (12) Vergani et al. (2011) (13) Page et al.(2007) (14) Yuan et al.(2010); (15) 
Marshall et al.(2011) (16) Golenetskii et al. (2010) (17) Nardini et al.(2014) (18) Tanvir et al.(2010) (19) Morgan et al.(2014) 
(20) Elliott et al.(2014) (21) Kelly et al.(2013) (22) Bellm et al.(2008) (23) Krim et al.(2007) (24) Pandey et al.(2010) (25) Rau 
et al.(2010) (26) Filgas et al.(2012) (27) Barthelmy et al.(2010GCN11023) (28) Hartoog et al.(2013) (29) Cucchiara et al.(2011) 
(30) Elliott et al.(2013) (31) Guidorzi et al.(2014) (32) Kruehler et al.(2013) (33) Vestrand et al.(2014) 

a The time intervals, for the fit results of the optical afterglow for each GRB; and the break time for the optical break samples, 
in units of ks. 

^The decay slopes before and after break time, note that although we also fit the rise slope for the GRBs which have the 
afterglow onset behaviors (such as GRB 071010A, 100901A), we did not summarize the information of the rise slope and peak 
time in our tables, since, in this paper, we focus on the decay behaviors. 

c Optical fiuence, which is calculated by integrating to flux over the time interval, in units of 10 — erg ern"^; the luminosity 
at the break time, in units of 10^ erg s — ^; and the isotropic energy release in the R-band, in units of 10^® erg. 

^The spectral index of the optical afterglow and the redshift of gamma-ray bursts. 

















Table 3 


The Relation Parameters of X-ray and Optical Afterglows 


GRB 

T x ,b/T 0 , b a 

Aai b 

AcK 2 b 

4a, h 

Ao:o b 

A/3 x b 

L X,b/ L R° 

®X,b/ B R C 

Fx/F 0 c 


Group I 

050319 

0.15i0.05 

-0.06i0.04 

-0.39=h0.20 

1.07i0.20 

1.40i0.00 

0.25i0.24 

144i30 

23i9 

42i7 

050408 

11.35i5.79 

0.37i0.03 

0.49=h0.20 

0.79i0.21 

0.67i0.14 


13i6 

237il02 

146i68 

050730 

1.15i0.15 

0.25i0.12 

1.24=1=0.05 

1.97i0.12 

0.98i0.03 

0.17i0.08 

164i29 

74il8 

217i28 

050922C 

0.51i0.26 

0.51i0.05 

-0.OliO.05 

0.53i0.06 

1.04i0.04 

0.02i0.16 

30il8 

21il4 

13i28 

060210 

27.00i37.64 

0.72i0.41 

0.19=1=0.14 

0.58i0.09 

1.12i0.42 

0.47i0.11 

17i27 

704ill95 

324i318 

060526 

0.78i0.18 

-0.10i0.60 

0.64=1=0.48 

2.26i0.77 

1.51i0.03 

-0.05i0.50 

6i2 

3i282 

3il 

060605 

0.35i0.05 

-0.48i0.06 

-0.40=1=0.15 

1.73i0.09 

1.65i0.13 

0.16i0.20 

35i9 

10i3 

7il 

060614 

0.76i0.05 

-0.04i0.05 

-0.03=1=0.05 

1.85i0.06 

1.84i0.03 

0.13i0.18 

28i2 


22i3 

060708 

12.00i6.31 

0.53i0.11 

0.46=1=0.05 

0.72i0.07 

0.79i0.10 

-0.OliO.28 

lil 

7i4 

7ill 

060714 

0.59i0.21 

0.23i0.11 

0.14=1=0.06 

0.88i0.12 

0.97i0.05 

0.19i0.23 

43il6 

36il6 

30il0 

060729 

1.23i0.18 

-0.OliO.10 

0.03=1=0.06 

1.27i0.02 

1.23i0.12 

-0.OliO.06 

7il 

10i2 

6il 

060912A 

2.18i2.15 

-0.14±0.34 

-0.19=1=0.05 

0.49i0.34 

0.54i0.02 

0.38i0.48 

3i3 

14il5 

18il6 

061021 

0.19i0.06 

-0.05i0.03 

-0.85=1=0.24 

0.55i0.04 

1.35i0.24 

-0.05i0.12 

113i36 

29il2 

33il6 

070411 

0.15i0.15 

0.00=h0.09 

-0.78=1=0.22 

0.54i0.24 

1.32i0.00 

-0.05i0.39 

105ill9 

15il0 

llil2 

070518 

4.83il5.07 

-0.01=h0.18 

-0.56=1=0.64 

0.65i0.64 

1.20i0.17 

0.52il.32 

5i8 

36il04 

llil7 

080310 

4.99i0.88 

0.13=h0.09 

0.51=1=0.06 

1.50i0.11 

1.12i0.03 

-0.15i0.23 

4il 

lli2 

7i3 

080413B 

0.27i0.05 

0.39=h0.01 

-1.31=1=0.12 

0.86i0.12 

2.56i0.01 

-0.07i0.30 

34i8 

18i6 

21i4 

080710 

1.74i0.51 

0.36=1=0.15 

0.39=1=0.14 

l.lli0.21 

1.08i0.01 

0.37i0.24 

3il 

5i2 

5i2 

080721 

0.16i0.01 

-0.20=1=0.01 

0.17=1=0.01 

0.85i0.02 

0.47i0.01 

0.15i0.06 

1794il34 

151il2 

180il0 

081029 

1.08i0.25 

-0.05=1=0.20 

0.62=1=0.17 

2.27i0.16 

1.60i0.21 

0.12i0.23 

4i 1 

5i2 

4il 

090618 

0.28i0.03 

0.02=1=0.01 

-0.00=1=0.06 

0.79i0.02 

0.81i0.06 

-0.07i0.06 

194i20 

59i8 

61i5 

090727 

112.67il48.33 

0.44=1=0.34 

0.82=1=0.30 

1.09i0.27 

0.71i0.37 

-0.17i0.61 



22ill 

091018 

0.OliO.00 

-0.59=1=0.06 

-0.43=1=0.02 

0.93i0.06 

0.77i0.00 

0.20i0.17 

2188i359 

14i3 

13i5 

091127 

1.07i0.55 

0.57i0.10 

0.10=1=0.13 

0.56i0.05 

1.03i0.16 

-0.03i0.14 

34il9 

62i47 

103i40 

Group II 

050401 



0.72=1=0.07 

0.92i0.07 


0.OliO.26 


749i472 

2184i888 

050502A 



0.56=1=0.05 

1.00i0.08 




0i9 

60il2 

050824 



0.38=1=0.13 

0.63i0.19 


0.63i0.42 


4i2 

8i4 

051221A 



0.48=1=0.07 

1.25i0.17 


0.10i0.27 


40i25 

lli4 

060111B 



0.46=1=0.16 

0.72i0.16 


0.04i0.38 



85i63 

060418 



0.30=1=0.06 

0.70i0.14 


-0.03i0.31 


10i5 

8i4 

060906 



0.98=1=0.23 

1.53i0.18 


-0.30i0.33 


143i51 

24i40 

060908 



0.47=1=0.07 

1.12i0.12 


-0.07i0.36 


32ill 

26i8 

060927 



0.82=1=0.46 

1.31i0.47 


0.10i0.43 


74i79 

8i7 

061121 



0.41=1=0.23 

1.07i0.03 


-0.13i0.13 


158i33 

37i25 

070318 



1.27=1=0.51 

1.56i0.51 


-0.20i0.50 


15i5 

6i7 

070420 



0.57=1=0.09 

1.21i0.06 


-0.15i0.19 



43i27 

071031 



0.91=1=0.48 

0.78i0.49 


-0.47il.02 


8il 

7il 

080810 



0.67=1=0.11 

0.70i0.18 


0.14i0.25 


4i3 

10il3 

081008 



0.89=1=0.14 

0.98i0.15 


-0.22i0.35 


61i22 

64i94 

081126 



1.17=1=0.06 

0.96i0.07 


0.39i0.25 



6i6 

081203A 



0.68=1=0.12 

0.98i0.12 


0.06i0.37 


7i2 

5i2 

090510 



1.17=1=0.15 

1.53i0.10 


0.38i0.22 


1540i894 

736i276 

100418A 



0.29=1=0.17 

1.68i0.11 


-0.05i0.34 


31i20 

17i67 

100621A 



1.20=1=0.10 

0.80i0.09 


0.72i0.33 



496i466 

100814A 



0.19=1=0.08 

1.59i0.06 


0.18i0.14 



261i230 

100906A 



1.05=1=0.80 

1.36i0.80 


0.09i0.18 


52il2 

18il4 

120119A 



0.67=1=0.18 

1.03i0.18 


0.31i0.40 


126i49 

61il9 

121217A 



0.66=1=0.05 

0.96i0.10 


-0.OliO.23 



55i96 

Group III 

050416A 



-0.43=1=0.32 


0.93i0.35 



631i388 

95i30 

050525A 



0.OliO.08 


0.75i0.07 



16i9 

54i22 

050801 



-0.13i0.06 


1.16i0.02 



5i2 

3i4 

051109A 



0.15i0.04 


0.65i0.09 



31i23 

56i26 

051111 



-0.49i0.70 


1.29i0.70 



2i2 

9i7 

061126 



-0.05i0.07 


0.67i0.09 



133i66 

166i54 

080413A 



-0.62i0.30 


1.18i0.30 



6i3 

15ill 

090426 



-0.21i0.07 


0.91i0.15 



55i29 

25i45 

130702A 



-0.06i0.29 


0.78i0.37 



33i34 

163il28 

Group IV 

050603 



-0.09i0.06 





5i3 

3i2 

050721 



0.00i0.03 






32i7 

050820A 



0.33i0.01 





24i3 

56i41 

060124 



0.45i0.02 





20i8 

91il8 

060206 



-0.07i0.02 





15i5 

4il 

061007 



0.06i0.01 





18il 

8i3 































Table 3—Continued 


GRB 

T x , b /T 0<b * 

Aai b 

AcK 2 b 

Aa x b 

Aa 0 b 

A/3 x b 

L X,b/ L R° 

S X ,b /Er c 

Fx/Fo° 

070125 



-0.01±0.19 





9±45 

3±3 

070311 



0.36±0.07 






9±6 

071010A 



-0.03±0.14 





8=1=1 

8=1=6 

071025 



0.51±0.04 





87=1=15 

22=1=32 

071112C 



0.45±0.02 





80=1=18 

42±56 

080319B 



0.31±0.00 





192±10 

76±2 

080804 



0.18±0.02 





107±35 

39=1=7 

080928 



-0.53±0.08 





11±11 

2=1=6 

081109A 



0.29±0.03 





55=1=18 

38=1=28 

090102 



0.28±0.01 





223=1=46 

175=1=19 

090323 



-0.15±0.21 





15±328 

3=1=8 

090328 



0.14±0.14 





11=1=31 

11=1=18 

090902B 



0.40±0.07 





50±55 

124±103 

090926A 



-0.06±0.04 





44=1=25 

2=1=1 

091029 



0.68±0.05 





272±72 

4=1=15 

100728B 



0.50±0.07 





44=1=25 

74=1=26 

100901A 



0.10±0.03 





8±1 

34=1=94 

101024A 



0.35±0.09 






120±105 

110205A 



0.20±0.01 





15=1=2 

7=1=7 

110213A 



1.17±0.09 





11±4 

23=1=89 

110918A 



0.29±0.05 





70=1=57 

112=1=63 

120404A 



0.32±0.12 





12±4 

7=1=66 

120815A 



0.28±0.04 





48=1=19 

40=1=54 

130427A 



0.25=1=0.00 





83±2 

83±1 


a The temporal ratio of break times (Group I samples). 

b The various values of pre-break decay (a: x i — a: 0 i)and post-break decay(o ; X) 2 — « 0) 2 )i the value of decay slope 
of pre- and post break of X-ray (Group I and Group II) and optical break samples (Group I and Group III); and the 
various values of the spectral index for X-ray break samples (Group I and Group II). 

c The ratio of break luminosity, isotropic energy release and flux at 1 hour between X-ray and optical band. 
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Table 4. Test the Temporal Indices for Various External Shock Models 


Total a Chromatic?* 3 

Injection 0 

Standard 01 

Standard 01 

Injection 0 

Chromatic?* 3 

External® 

Standard® 

wind 

wind 

wind 

ISM 

ISM 

ISM 

ISM and wind 

ISM and wind 

< -0.5 

(-0.5, -0.25) 

(-0.25, 0.0) 

(0.0, 0.25) 

(0.25, 0.5) 

> 0.5 

(-0.5, 0.5) 

(-0.25, 0.25) 


Group Samples 


Group I 


Aa(a x ,l — «o,l) 

24 

1 

1 

8 

5 

4 

5 

18(75%)/23(96%)* 

13(54%)/16(67%) 

Aa(a x ,2 — a o,2) 

24 

3 

3 

5 

5 

3 

5 

16(67%)/21 (88%) 

10(42%)/16(67%) 

Group II 

Aa(a X| i — £* 0 , 2 ) 

24 

7 

11 

3 

2 

1 

0 

17(71%)/20(83%) 

5 (21%)/10(42%) 

A «(« x ,2 — £* 0 , 2 ) 

24 

0 

0 

0 

1 

7 

16 

8(33%)/14(58%) 

1(4%)/7(29%) 

Group III 

A«(a X| 2 — £*o,l) 

9 

0 

0 

0 

0 

3 

6 

3(33%)/6(67%) 

0(0)/l(ll%) 

Aa (« x ,2 — £* 0 , 2 ) 

9 

1 

2 

4 

2 

0 

0 

8(89%)/9(100%) 

6(67%)/8(89%) 

Group IV 

A «(« x ,2 — £* 0 , 2 ) 

30 

1 

0 

6 

7 

14 

2 

27(90%)/29(97%) 

13(43%)/20(67%) 

Epoch Samples 

epoch 1 

Aa(a x — ot Q ) 

50 

6 

7 

9 

7 

10 

11 

33(66%)/42(84%) 

16(32%)/23(46%) 

epoch 2 

Aa(a x — ct Q ) 

97 

8 

13 

15 

15 

22 

24 

65(67%)/76(78%) 

30(31%)/45(46%) 

epoch 3 

Aa(a x — ct Q ) 

80 

12 

6 

13 

13 

15 

21 

47(59%)/60(75%) 

26(33%)/37(46%) 

epoch 4 

Aa(a x — a Q ) 

35 

4 

2 

5 

7 

12 

5 

26(74%)/28(80%) 

12 (34%)/19(54%) 


a Total numbers of data points for each Group or epoch samples. 

^Numbers of the data points for which are Aa < —0.5 or Act > —0.5, which do not fall into our external shock model ranges. 
c Numbers of the data points for energy injection model, ISM and wind-like medium. 

^Numbers of the data points for standard afterglow model, ISM and wind-like medium. 

e Numbers of the data points (brackets correspond with the fraction) for satisfied with the external shock and standard afterglow model ranges for each Groups 
or epochs samples. 

^Not consider error bar/consider error bar. 


Table 5. The Probability of the K-S Test for Our Epochs Samples 


epochs 

Ck — S ( a x,obs i a x,Sin 

o a 

^K-s(“o,obs. “o,Sin 

u b 

epoch 1 
epoch 2 
epoch 3 
epoch 4 

0.02 

8.36x 10 —3 
0.24 

0.03 


1.32 x 10 —10 
2.28xl0~ 9 
7.32x 10 —7 
0.49 



a The probability of the K-S test for observed decay slope and simu¬ 
lated decay slope for the X-ray band. 

b The probability of the K-S test for observed decay slope and simu¬ 
lated decay slope for the optical band. 


























Table 6. The Closure Relation of Afterglows in Gamma-ray Bursts 


oc(q = 1) 
Non-Injection 


a(0 < q < 1) 
Injection 


«(/3)(« = 0 

Non-Injection 


c(/3)(0 < q < 1) 
Injection 


a(P) 

Jet 


ISM, Slow cooling 


. (2+q)P 

(2+q)P 

2 _ 


v < v m 

n < V ^ v c 
v > v c 


3(1 -p) 


2-3 p 


8 —5g 

(6-2p)-(p+3)g 

(4-2p)-(p + 2)g 

_4_ 


' —P 
' ~P 


_ 3/3-1 

LZ_2_ 


« = (<?- 1) + 

a = (q- - 1) + 

— g~ 2 

__ 2 _ 


a: = 2/3 + 1 
oc = 2/3 


ISM, Fast cooling 


I 

r 2 
f (2 — q)/3 

(2+g^/3 

_ 2 _ 


v <v c 
2 < v < v n 
v > Vm 


2-3 p 


8 —7g 
6 

2 —3g 

(4-2p)-(p + 2)g 

_4_ 


“ — 2 

a = £ 

2 

_ 3/3-1 

1_Z_2_ 


: (g - 
: (g- 


, - q~ 2 

__ 2 _ 


Wind, Slow cooling 


- TJ=qT - 

(2-2p)-(p + l)g 

(4-2p)-(p + 2)g 

_4_ 




(2+ 2 g)/3 
(2 + g)/3 

I_2_ 


v < 

rn < v < v c 
v > v c 


0 

1 — 3p 
4 

2-3 p 


' ~P 
' -P 


. _ 3/3 + 1 
— 2 
_ 3/3-1 

LZ_2_ 


OL = -k - 


. - q-2 

__ 2 __ 


a: = 2/3 + 1 
a = 2/3 


Wind, Fast cooling 


- EH - 

3 

2 —3g 

(4-2p)-(p + 2)g 

_4_ 


_ 1-/3 
* ~ 2 
— !-/3 
* ~ 2 
_ 3/3-1 

LZ_2_ 


_ q (2-qJB- 
, _ 1 _ (2-q)/3 


v <v c 
v c < ^ 

^ > fm 


2 —3p 


a = f — ■ 

g-2 , (2 + g)/3 

__ 2 _ 1 _ 2 _ 


> = 2.5(p = 2,p = 3), g = 0 


ISM, Slow cooling 


i' < vm 0.5(0.5, 0.5) 1.3(1.3, 1.3) 0.5(0.5, 0.5) 0.3(0.3, 0.3) 

rn<u<v c -l.l(-0.75, -1.5) 0.25(0.5,0.0) -2.5(-2.0, -3.0) -0.75(-0.75, -0.75) 

v > v c -1.4(-1.0,-1.75) -0.25(0.0, -0.5) -2.5(-2.0,-3.0) -1.25(-1.25,-1.25) 


ISM, Fast cooling 


v<v c 0.2(0.2, 0.2) 1.3(1.3, 1.3) 

z<v<vm -0.25(-0.25,-0.25) 0.5(0.5,0.5) 

v > vm _—1.4( —1.0, -1.75)_-0.25(0.0, -0.5) 


0.3(0.3, 0.3) 

—0.5(—0.5, -0.5) 

—1.25( —1.25, -1.25) 


Wind, Slow cooling 


v < v m 0.0(0.0, 0.0) 0.3(0.3, 0.3) 0.5(0.5, 0.5) 0.3(0.3, 0.3) 

rn<u<v c -1.6(-1.25,-2.0) -0.75(-0.5,-1.0) -2.5(-2.0,-3.0) -0.75(-0.75,-0.75) 

v > v c _—1.4( —1.0, -1.75)_-0.25(0.0, -0.5) -2.5(-2.0,-3.0) -1.25(-1.25,-1.25) 


Wind, Fast cooling 


v < v c 
v c < v < Vm 
v > Vm 


—0.7(—0.7, -0.7) 

—0.25(—0.25, -0.25) 
—1.4( —1.0, -1.75) 


—0.3(—0.3, -0.3) 
—0.5(—0.2, -0.5) 
-0.25(0.0, -0.5) 


0.3(0.3, 0.3) 

—0.5(—0.5, -0.5) 

—1.25( —1.25, -1.25) 


Notes: This table provides the ranges of relations in the temporal index and the spectral index that are expected from synchrotron emission with or without energy injection, and the post 
jet-break case (Zhang Sz Meszaros 2004; Zhang et al. 2006; Sari et al 1999; Panaitescu et al. 2006). Here p is the electron index and q is the luminosity index. When q = 1 non-injection is 
the case and when q = 0 energy injection is the case; and calculated numerical values for each cases with p = 2.5 p = 2.0, p = 3.0, q = 1 and q = 0. 
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Table 7. Diagnosis of the Properties of Afterglow Break Behaviors. 


GRB 

“1 < 

t*2 > l.s a 

Closure Relation?' 3 

Spectral Evolution 0 

Spectral Break 0 

Injection Break? 0 

Jet break? 0 


Group I 

050319 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

050408 

No/Yes 

Yes/No 

Yes/Yes 

No 

No 

No/Yes 

Yes/No 

050730 

Yes/Yes 

Yes/No 

No/Yes 

Yes 

No 

Yes/Yes 

Yes/No 

050922C 

Yes/Yes 

No/No 

Yes/Yes 

No 

No 

Yes/Yes 

No/No 

060210 

No/Yes 

No/No 

Yes/Yes 

No 

No 

Yes/Yes 

No/No 

060526 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

060605 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

060614 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

060708 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

060714 

Yes/Yes 

No/No 

Yes/Yes 

No 

No 

Yes/Yes 

No/No 

060729 

Yes/Yes 

No/No 

Yes/Yes 

No 

No 

Yes/Yes 

No/No 

060912A 

Yes/No 

No/No 

Yes/Yes 

No 

No 

Yes/Yes 

No/No 

061021 

Yes/... 

No/Yes 

Yes/Yes 

No 

No 

Yes/... 

No/Yes 

070411 

Yes/... 

No/Yes 

Yes/Yes 

No 

No 

Yes/... 

No/Yes 

070518 

Yes/Yes 

No/No 

No/Yes 

No 

No 

Yes/Yes 

No/No 

080310 

Yes/Yes 

Yes/No 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/No 

080413B 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

080710 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

080721 

Yes/No 

Yes/No 

Yes/Yes 

Yes 

No 

Yes/No 

Yes/Yes 

081029 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

090618 

Yes/Yes 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/Yes 

090727 

Yes/Yes 

Yes/No 

Yes/Yes 

No 

No 

Yes/Yes 

Yes/No 

091018 

Yes/No 

No/Yes 

Yes/Yes 

No 

No 

Yes/No 

No/Yes 

091127 

No/Yes 

Yes/No 

Yes/Yes 

No 

No 

No/Yes 

Yes/No 

Group II 

050401 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

050502A 

Yes/... 

Yes/No 

Yes/Yes 


No 

Yes/... 

Yes/No 

050824 

Yes/... 

No/No 

No/Yes 

Yes 

No 

Yes/... 

No/No 

051221A 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

060111B 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

060418 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

060906 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

060908 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

060927 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

061121 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

070318 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

070420 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

071031 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

080810 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

081008 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

081126 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

081203A 

No/... 

Yes/Yes 

Yes/Yes 

No 

No 

No/... 

Yes/Yes 

090510 

Yes/... 

Yes/No 

Yes/Yes 

Yes 

No 

Yes/... 

Yes/No 

100418A 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

100621A 

No/... 

Yes/No 

No/Yes 

Yes 

No 

No/... 

Yes/No 

100814A 

Yes/... 

Yes/Yes 

Yes/Yes 

No 

No 

Yes/... 

Yes/Yes 

100906A 

Yes/... 

Yes/No 

Yes/Yes 

No 

No 

Yes/... 

Yes/No 

120119A 

No/... 

Yes/No 

Yes/Yes 

No 

No 

No/... 

Yes/No 

121217A 

Yes/... 

No/No 

Yes/Yes 

No 

No 

Yes/... 

No/No 

Group III 

050416A 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 

050525A 

.../Yes 

Yes/Yes 

Yes/Yes 



.../Yes 

Yes/Yes 

050801 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 

051109A 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 

051111 

.../Yes 

Yes/Yes 

Yes/Yes 



.../Yes 

Yes/Yes 

061126 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 

080413A 

.../Yes 

No/Yes 

Yes/Yes 



.../Yes 

No/Yes 

090426 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 

130702A 

.../Yes 

No/No 

Yes/Yes 



.../Yes 

No/No 


a X-ray/optical. We defined the shallow decay by the ’closure relation’ with different spectral regimes, 3 GRBs lack the optical spectral 
information; we choose the jet break with the criterion of a 2 > 1.5. 

^Test the data with the ’closure relation’ (ck 2 — ^ 2 )- 

c The spectral breaks required Acc = 0.25 and A/3 = 0.5; only for X-ray break samples, since no time-resolved spectral analysis for optical 

























band available; and statistic break typi 
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Fig. 1.— The cartoons showing onr four samples which, with different shapes with or without 
breaks, comparing X-ray and optical afterglow light-curves and the models used. Note that 
Group II GRBs do not always point towards a chromatic behavior. This is because in many 
GRBs there is no optical data around the X-ray break time, so that the existence of an X-ray 
break around the same time is not ruled out. Again some GRBs in Group III do not have 
X-ray data around the optical break time, so that they are not necessarily chromatic. We 
therefore plot the post-break segment for Group II and the pre-break segment for Group III 
by the dot-dashed lines. 






Fig. 2.— Examples of the best afterglow light-curve fittings, that distinguish four shapes 
with and without breaks for our samples. 
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Time since trigger (s) 








Time since trigger (s) 





Fig. 3.— The cases in which breaks were detected in their afterglow light-curves both in the 
X-ray and the optical band. 
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Fig. [3]— Continued 
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Fig. 4.— The cases in which a breaks was detected only in the X-ray band but not in the 
optical band in their afterglow light-curve. 
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Time since trigger (s) 









Fig. 0]— Continued 
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Fig. 5.— The cases in which a breaks was detected only in the optical band but not in the 
X-ray band in their afterglow light-curves. 
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Fig. 6.— The cases that did not have detected break behavior in neither the X-ray band nor 
optical band in their afterglow light-curves. 
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Time since trigger (s) 







Fig. IS]— Continued 
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Time since trigger (s) 





Fig. |6]— Continued 
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Fig. 7.— The distributions of properties for our cases comparing the X-ray with optical 
emission, including the temporal index, the break time and break luminosity (for the break 
cases), and flnence and isotropic energy release. 




WJ 


Fig. 8.— Pair correlations comparing the X-ray with optical emission for Group 1 cases, 
including the change of X-ray and optical pre- and post-break decay slopes and break time. 
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Fig. 9.— The time interval of the optical detection as a function of the time interval of the 
X-ray detection for all afterglow cases. The dashed lines are the equal lines. 
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a x 



a x 




a x a x 


Fig. 10.— The power-law indices for the light-curve decay. The panel from the top to 
bottom show Group I, II, III and IV,respectively. The black solid line indicates where are 
the X-ray and optical temporal indices are equal and the red and green solid line indicates 
where a x = a Q ± 0.25, which is expected by standard forward shock model with ISM and 
wind-like environment. The red and green dashed lines are a x = a Q ± 0.5, which is expected 
by the energy injection case for a ISM and a wind-like environment. 
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Fig. 11.— The power-law indices of spectra. The panel from the top to bottom show Group 
I, II, III and IV, respectively. The black solid line indicates where are the X-ray and optical 
spectral indices are equal (Spectral regime I: v m < v Q < u x < u c ) and the black dashed line 
indicate where /3 X = (3 0 + 0.5 (Spectral regime II: is m < u Q < u c < z/ x ). 
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Fig. 12. — The distribution of the temporal and spectral indices for both the X-ray and 
optical bands for our afterglow cases, (a): The lines are the best Gaussian fitting with the 
center value a x ,i( c )=0.68± 0.03, a o ,i(c)=0.44± 0.01 for the pre-break decay slopes and (b): 
«x, 2 ( c ) = 1-48± 0.02, a 0 , 2 (c) = l-23± 0.03 for the post-decay slopes, respectively. The center 
value of all cases are consistent with case of ISM medium and slow cooling, u m < u Q < 
v c < iv (c): The best Gaussian fitting for spectral are /S x ,2(c)=0-99± 0.01, /S o ( c )=0.69± 0.01, 
respectively. This is generally consistent with the model with ISM medium and slow cooling, 
< Vo < v c < J'x- 
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a x ,10 4 s-10 5 /(1+z)s a x ,>10 5 /(1+z)s 


Fig. 13.— As Figure [TO], the external shock model is compared to the power-law indices 
measured over different time intervals for the 4 groups. The panel from the left to right 
show < 10 3 s, 10 3 — 10 4 s, 10 4 — 10 5 s, > 10 5 s for each bin, respectively. Note that we have 
corrected them to the rest frame. The black solid points are observe data and red hollow 
points are simulate data. 
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Fig. 14.— The temporal decay index a versus spectral index f3 along with the ’closure 
relations’ of the external forward shock models for all our cases. Regime I {y > max(v m , z/ c )), 
Regime II (z/ m < v < u c ). The thick solid lines and solid shaded regions are the closure 
relation of the pre- and post-break in spectral regime I. The lower and upper boundaries of 
the regions are defined by the ’closure relation’, respectively. Similarly, the thick dashed lines 
and hatched regions are the spectral regime II. The solid round and triangle dots of black 
and blue are the temporal decay index of pre- and post-break for our total break samples 
(X-ray: Group I+Group II; Optical: Group I+Group III.); the half-solid diamond of pink is 
the temporal index for our non-break samples (X-ray: Group III+Group IV; Optical: Group 
II+Group IV.). (a) The case of ISM model for XRT data, (b) The case of wind model for 
XRT data, (c) The case of ISM model for optical data, (d) The case of wind model for 

nntirfll rla.t.p. 
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Fig. 15.— The temporal decay indices of pre-break oi\ versus post-break a 2 for our break 
cases, obtained the optical break samples (Group I+Group III) and the X-ray break samples 
(Group I+Group II). The three solid lines are the various external shock models which 
invoking energy injection with index q= 0. 
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q 




Fig. 16. — The relation of Aa, q and p indices. The distributions of Aa = a 2 —a i, luminosity 
injection index q and electron spectral index p (Aa calculated based on 44 X-ray break cases 
(Group 1+ II) and 33 optical break cases (Group I+III), and the q index is derived from 
the shallow decay segments, and the best Gaussian function fitting with the center value 
g x ( c )=0-46± 0.03 and g o ( c )=0.13± 0.11, respectively. The p index is derived from the post¬ 
break decay segments for all the cases, and the best Gaussian fitting with a center value 
Px(c)=2.63±0.02 and p o ( c )=2.58±0.04 for X-ray and optical. The solid and dashed lines are 
the best Gaussian fitting for X-ray and optical, respectively. The right panel is Aa — p, q 
planes, in panel (c), the black, red and blue line are the q— 0, g=-0.5 and <7=0.5, respectively. 
The black and red points are the X-ray and optical data, respectively. In panel (d), the black 
solid line and red dash line are p —2 and p —3 for different external shock models. 



























Fig. 17. — The cooling frequency crossing the X-ray and optical band. The change in 
temporal slope A« = «2 — ot\ for both the X-ray (Group I+II) and optical (Group 1+ III) 
break cases and the spectral index A/3 X = /3 X>2 — /3 Xi i for the X-ray break cases (Group 
I+Group II), and the best Gaussian function fitting with the center value /3 x ,i(c)—0-95± 0.03 
and /3 X]2 (c)=0.98± 0.01, respectively. In panel (c) shows the relation of spectral slopes of pre 
and post break segments for X-ray break cases, the black solid line indicates where are the pre 
and post spectral indices are equal and the black dashed line indicate where /3 X) 2 = /3 x ,i ±0.5 





























